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BASIC REQUIREMENTS FOR MANUSCRIPTS 


This Journal represents effort the Society deliver information the 
reader with the greatest possible speed. this end the material herein has 
none the usual editing required more formal publications. 


Original papers and discussions current papers should submitted the 
Manager Technical Publications, ASCE. Authors should indicate the technical 
division which the paper should referred. The final date which dis- 
cussion should reach the Society given footnote with each paper. Those 
who are planning submit material will expedite the review and publication 
procedures complying with the following basic requirements: 


Titles should have length not exceeding characters and spaces. 
summary should accompany the paper. 


The manuscript ribbon copy and two copies) should double-spaced 
one side paper. Papers that were originally prepared for 
oral presentation must rewritten into the third person before being submitted. 


full name, Society membership grade, and footnote reference 
stating present employment should appear the first page the paper. 


Mathematics are reproduced directly from the copy that submitted. 
Because this, necessary that capital letters drawn, black ink, 
high (with all other symbols and characters the proportions dictated 
standard drafting practice) and that line mathematics longer than 
Ribbon copies typed equations may used but they will proportionately 
smaller the printed version. 


Tables should typed (ribbon copies) one side 11-in. 
paper with invisible frame. Small tables should grouped 
within this frame. Specific reference and explanation should made the text 
for each table. 


paper within invisible frame that measures the caption 
should also included within the frame. Because illustrations will reduced 
69% the original size, the capital letters should high. Photographs 
should submitted glossy prints size that less than 
Explanations and descriptions should made within the text for each illustration. 


Papers should average about 12,000 words length and should 
longer than 18,000 words. approximation, each full page typed text, 
table, illustration the equivalent 300 words. 


Further information concerning the preparation technical papers con- 
tained the “Technical Publications Handbook” which can obtained from 
the Society. 


Reprints from this Journal may made condition that the full title 
the paper, name author, page reference (or paper number), and date 
publication the Society are given. The Society not responsible for any 
statement made opinion expressed its publications. 


This Journal published monthly the American Society Civil 
Engineers. Publication office 2500 South State Street, Ann Arbor, Michigan. 
Editorial and General Offices are West Street, New York 18, New York. 
$4.00 dues are applied subscription this Journal. Second- 
class postage paid Ann Arbor, Michigan. 


Subject and author indexes, with abstracts, are published the end each 
year for the Proceedings ASCE. The index for 1958 was published Proc. 
Paper 1891; indexes for previous years are also available. 
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MATRIX SOLUTION BEAMS WITH VARIABLE MOMENTS 
INERTIA 


SYNOPSIS 


Using finite differences and matrix algebra, new formulas are developed 
which permit the direct computation bending moments, moment influence 
coefficients, deflections, and deflection influence coefficients simply sup- 
ported beams with variable moments inertia. These formulas are explicit 
the desired quantities and are especially suitable for use digital com- 
puter programs. 


Matrix Formula for Bending Moment 


Consider simply supported beam subjected arbitrary transverse 
loading function q(x), shown Fig. The bending moment the 
beam, governed the differential equation, 


(1) 
together with the two boundary conditions, 


Let divide the span length the beam into equal intervals each; 
and designate the equally spaced points Points n-2, n-1, and 
Points and are the left and right support respectively. Eq. (1) may 
expanded into finite difference equations. 


Note: Discussion open until March 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2218 
part the copyrighted Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. October, 1959. 
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FIG 


One such equation may written for each the points 2,..., n-1, giving 
total n-1 equations. 


M; 2M, My 4) i= 1253900 (3) 


The boundary equations become, 


After substitution Eq. (4) into Eq. (3) order eliminate and 
have n-1 simultaneous linear algebraic equations containing n-1 unknown 
moments Mj. Writing out these equations full, have, 
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(5) 


The set n-1 simultaneous equations can condensed into single matrix 


equation. 


where [B] square matrix and and are two column matrices, de- 


fined the following: 


(6) 
(7) 
(8) 
(9) 
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The basic differential equation (1) together with the boundary conditions (2) 
are now replaced single matrix equation (6). solve Eq. (6) for the 


moment matrix pre-multiply both side the equation the inverse 
matrix 


the matrix [B]~1 appears frequently the following paragraphs, 
shall designate matrix 


The matrix independent the load matrix and can evaluated once 


and for all. fact, easy show matrix algebra that matrix 
given the following: 


(12) 
trical 


particular, matrix for three different values are given below. 


For 


For 


For 10, 


Eq. (11) the general solution the bending moment simple beam 
matrix form. Given any loading condition, the multiplication the two known 
matrices and gives once the values the bending moments the n-1 
equally spaced points along the beam. Eq. (11) not new method 
learned and mastered, but formula explicit the unknown bending 
moments. 

When the load function series concentrated loads, represented 
exactly the load matrix this case, Eq. (11) gives the exact solution 
moments the beam. When the load function cannot represented exactly 
the load matrix i.e., the case nonuniformly distributed load, 
the bending moments obtained Eq. (11) are only approximate. The approxi- 


mation, course, can improved without limit choosing larger values 


Example 


Compute the bending moments simple beam due the set loads 
shown Fig. 


The element the load matrix qh, represents the magnitude the con- 
centrated load the interval 


From Eq. (11), have, 
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206 

228 
Example 


Compute the bending moments simple beam due uniformly dis- 
tributed load applied over the entire beam. this case the load matrix 
simply the following. 


The bending moments are computed Eq. (11). 


125 


can easily verified that the elements the moment matrix are the 
exact bending moments the nine points. 


Moment Influence Coefficient Matrix 


The moment influence coefficient, defined the bending moment 


ai 
105 
105 


the beam Point due unit load applied Point the beam divided 
into equal segments, there will (n-1)2 these coefficients. Arranging 
them square array, have the moment influence coefficient matrix. 
compute these coefficient, apply unit load Points succes- 


sively. Carrying out the computation one operation, use unit matrix 
the load matrix. 


Eq. (11), compute the moment influence coefficient matrix. 


Thus seem that the moment influence coefficient matrix simply equal 
matrix For example, the third row matrix 


gives the bending moments Point due unit load applied Points 
successively. 


Matrix symmetrical matrix. That is, 

(14) 

This direct consequence the Maxwell’s Reciprocal Theorem which 

states that “the moment due unit load equal the moment 

The moment influence coefficient matrix simple beam may useful 
the analysis truss members under live which case the panel 
length the truss should taken 


Matrix Formula for Deflections 


The deflection, simple beam with variable moments inertia 
governed (see Fig. 1), 


(15) 


Again expanding Eq. (15) into finite difference equations and taking into 
consideration Eq. (16), have, 
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The n-1 equations Eq. (17) are condensed into one matrix equation. 


i=l, n-l 


where matrix has been defined previously Eq. (7). The deflection 


matrix and the Mh/EI] are both column matrices. 


The elements the matrix Mh/EI] represent the area the familiar M/EI 
diagram within the interval Solving Eq. (18) for the deflections, have, 


[4] 


where matrix the same matrix previously defined Eq. (12). With 
(21), one can compute the deflections direct substitution. 


Example 


Compute the deflections simple beam with constant moment inertia 
subjected two equal and opposite moments applied its ends. 


(19) 
(20) 


18) 


19) 


(20) 
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Example 


Calculate the deflections nine points along the beam carrying the loads 
shown Ex. The flexural rigidity the beam given Fig. 
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206 103 
228 
232 116 
130 
282 
8,154 
13,328 
9,776 
5,328 


The matrix multiplications required this example can easily carried out 
desk calculator. 


IV. Deflection Directly From Load 


some instances, desirable compute the deflections beam with 
variable moments inertia directly from the given load matrix Utilizing 
previous results, the new matrix formula derived. 


(EI y")" 


where the matrices and have been defined Eqs. (19), (12) and (9). 
The flexural rigidity matrix [h/EI] diagonal matrix. diagonal matrix 
square matrix having zeroes for all elements except for the diagonal 
elements. 


| 
| 
tg eee 


ith 
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23) 


(23) the most important equation derived this paper. With it, one 
can calculate the deflections from given loads and moments inertia 
direct matrix multiplications. 


Example 


Calculate the deflections the beam shown Ex. means Eq. (22). 
The given data are the load matrix and the flexural rigidity matrix 


Substituting the above matrices into Eq. (22), 


(22) involves three matrix multiplications. Because matrix multiplica- 
tions are associative, the order which these multiplications are carried out 
immaterial. 


Deflection Influence Coefficient Matrix 


The deflection influence coefficient, defined the deflection the 
beam Point due unit load applied Point 


(24) 


Written matrix notation, Eq. (24) becomes, 


where and are the deflection and load matrix respectively, and the 
deflection influence coefficient matrix. The elements matrix are 
deflection influence coefficients. (The deflection influence coefficients are 


also known the flexibility influence coefficients.(6)) Comparing Eq. (25) 
with Eq. (22), seen once that, 
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Eq. (26) the general formula deflection influence coefficient matrix for 
beams with variable moments inertia. For beams with constant moment 
inertia, Eq. (26) may simplified. 


[1] 


Several [A]2 are given Appendix for easy reference. 


Example 


Find the deflection influence coefficient matrix for beam constant 
moment inertia. Using the value given Appendix for 10, 


285 480 595 640 625 560 320 165 
480 880 1120 1220 1200 1080 880 620 320 

595 1120 1505 1680 1625 1520 1245 880 455 

640 1220 1680 1960 2000 1840 1520 1080 560 

[c] 625 1200 1675 2000 2125 2000 1675 1200 625 
560 1080 1520 1840 2000 1960 1680 1220 640 

455 880 1245 1520 1675 1680 1505 1120 595 

320 620 880 1080 1200 1220 1120 880 

455 560 625 640 595 480 285 


For example, the coefficient for the deflection the midspan due unit 
load the mid-span is, 


Compared with the exact answer the error about 2%. 


Example 


Determine the deflection influence coefficient matrix for the beam shown 
Fig. 


FIG 


For the sake simplicity, shall divide the beam into five equal segments 
and carry out the computations detail. Eq. (26), 
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Thus with little effort, mental computational, have obtained sixteen de- 
flection influence coefficients. The results are, course, not very precise, 
only five segments were taken. for example, has error about 
10.5%. 


CONC LUSIONS 


There are many excellent methods for the calculation beam deflections. 
Among the most well-known are the ate beam method, (7) Newmark’s 
method, (8) and the relaxation method.(L With any these methods, one 
can compute beam deflections systematically and efficiently. Their only 
drawback, can called drawback, that one must learn and memorize 
the particular procedure involved. Unfortunately the procedures can easily 
forgotten, unless one uses them fairly regularly. This probably the main 
reason why Newmark’s method did not have the popularity rightfully de- 
serves among the engineers. 

With the matrix formulas developed this paper, the deflection computa- 
tions are reduced mere routine. All one needs know are the formulas 
and the rule for matrix multiplication. would possible, for example, 
explain how use the formulas engineering aide in, say, half hour. 
Thereafter, will able crank out deflections desk calculator 
fast anyone. 

None the existing methods are particularly suitable for use digital 
computer, because the long and involved programming required. With the 
matrix formulas developed this paper, one can write short, possibly the 
shortest, subroutine for beam deflections any digital computer. The 
matrix can fed into the computer part the subroutine. For medium 
capacity machines, such the 650, G-15, LGP-30, entirely feasible 
desirable let the computer generate matrix internally, for any value 

The extension the matrix formulas compute stiffness matrix, buckling 


load, and the natural frequency beams with variable moments inertia will 
treated subsequent paper. 
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116 


136 114 


144 171 


640 625 


1220 1200 
1680 1675 
1960 2000 
2000 2125 
2000 
1520 1675 


1080 1200 
455 560 625 640 


455 

880 
1520 
1675 
1680 
1505 
1120 


595 480 


ASD 


[4] 
136 
285 560 320 165 
480 1080 620 320 
595 1120 150 1520 880 
640 1220 168 1080 560 
560 1080 152 1960 1220 640 
455 880 1680 1120 595 
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REINFORCED CONCRETE FOLDED PLATE CONSTRUCTION 


Boyd Anderson,2 ASCE and 
Harold 


This paper one the group papers which were the bases for oral 
presentations the Joint ASCE-IABSE Meeting New York, October 1958. 


The entire group, including those published CIVIL ENGINEERING, will 
reprinted one volume. 


SYNOPSIS 


This paper presents design methods and discusses experience with the 
construction reinforced concrete folded plate thin shells for the roofs and 
floors hangars, auditoriums, warehouses and other buildings. Quantities 
materials and costs are compared with other types structures. 


INTRODUCTION 


Reinforced concrete folded plates, they are sometimes called, pris- 
matic shells hipped plates, provide useful and economical method con- 
struction for roof and floor systems wide variety structures. 
often competitive with other construction methods for short spans and has 
proven exceptionally economical where relatively large spans are needed 
for auditoriums, gymnasiums, industrial buildings, hangers, department 
stores and parking garages. 

The materials required are usually much less than needed for flat slab, 
beam and slab, other conventional systems and are little more than 


Note: Discussion open until March 1960. Separate Discussions should submitted 
for the individual papers this symposium. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2219 
part the copyrighted Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. October, 1959. 

Partner, Ammann Whitney, Cons. Engrs., New York, 

Partner, Ammann Whitney, Cons. Engrs., New York, 

Engr., Ammann Whitney, Cons. Engrs., New York, 


ll 
Kod 
3 
i 


i 


required for continuously curved shells, with the advantage utilizing rela- 
tively simple formwork. many cases, concrete folded plates will provide 
fire resisting structures substantially less cost than structural steel 
framing. 

Folded plate structures were first constructed Germany during the 
1920’s for use coal bunkers and similar structures where the ratio span 
width plate relatively small. Technical papers the subject began 
appearing about 1930. 

The first application America may have been warehouse San 
Francisco designed Nishkian, Consulting Engineer, and built about 
1935. This was about feet wide and feet long. The roof consisted 
three ribbed concrete flat segments with the center plate horizontal and the 
roof load carried the end walls the two sloping side plates. re- 
ported that number similar structures were subsequently built 
California. 

recent years, concrete folded plates have been used for large variety 
structures and are increasing importance the building industry the 
basis new system construction. 


Applications Folded Plates 


Folded plate floor roof construction consists series repeated 
units, each which formed two more flat plates intersecting 
angle. The plates act continuous slab transversely and beams their 
own planes longitudinally. 

The angle the sloping plates usually made 45° less with the hori- 
zontal order facilitate placing the concrete without top forms. The 
simplest form two-plate V-shaped unit. The disadvantage this form 
that may provide inadequate area concrete the top and bottom 
resist the compressive flexural force permit placement the reinforcing 
steel. 

The most useful and practical form unit developed adding hori- 
zontal plates between the sloping plates the top and bottom sufficient 
width keep the concrete compressive stress within proper limits and 
accommodate the necessary reinforcing steel. The majority the following 
examples are this type. 

Sometimes the bottom plate placed vertically, forming downward 
projecting beam instead horizontal flange. This has the disadvantage 
greater height structure and greater form cost. 

The use inclined plates which extend from top bottom the unit 
usually more economical both cost construction and design than more 
complicated forms with intermediate folds. 

Tapered triangular plates may also used where desired. the 
Illinois dome, (Fig. tapering makes possible the use folded plates the 
circular structure without excessive convergence and narrowing the folds. 
the roof the American Concrete Institute office building(1) they are used 
for aesthetic reasons. 

The folded plate units are usually symmetrical where the solid roof 
floor desired, but unsymmetrical units may used provide north-light 
roofs. 


se 
4 Ou 
cc 
S} 
ve 
a 
Ss. 
1 


pan 


the 


PLATE CONSTRUCTION 


number examples folded plate construction will illustrated 
the following pages and drawings are presented showing typical cross 
sections. 

The roofs the social hall and chapel the Temple Beth Sholem 
Miami (Fig. are formed four-inch concrete slabs corrugated height 
The spans were feet and feet respectively and the cost was 
comparable the cost much shorter spans using conventional framing. 

The roof the Sears Roebuck Co. store Tampa also the simple 
type corrugated construction (Fig. 1). consists two 125 foot continu- 
ous spans with cantilevers the ends. The strength and rigidity the 
construction indicated the fact that the roof supports steel and con- 
crete mezzanine floor with suspenders foot intervals leaving the 125 
foot spans clear the first story. 

Temple Mishkan Tefila, Newton, Massachusetts (Fig. 12), has foot 
span. The horizontal plate has been omitted the bottom forming units con- 
sisting V’s separated narrow horizontal plates the top. There 
vertical beam projecting below the forming the bottom 
flange. The total height the section feet. 

The floors the Wayne University conference building (Fig. are the 
simple type with level slab poured across the top forming the floor sur- 
face. The hollow triangular spaces were used for ventilating ducts. 

The folded plate roof the cold storage warehouse for the Central Ware- 
house Corporation Albany (Fig. 10) were also leveled with top slab 
the interior spans were used cold air ducts for refrigeration. The 
operation the plant has been very satisfactory. 

Another example the use folded plate construction with top levelling 
slab the central portion the TWA hangar Kansas City (Fig. 11). The 
upper floors were suspended from the roof with two foot spans, leaving the 
ground floor open for the full width 100 feet. The folded plate construction 
provided shallower and more economical framing than could obtained 
using beams and girders. 

Folded plate roofs have been used several high schools Milwaukee 
for gymnasiums, auditoriums, and other portions (Fig. 8). Examples 
these are the John Muir High School (Eschweiler Eschweiler, Architects), 
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ing 
Fig. Sears Roebuck Co. Store, Tampa, Florida 


the John Marshall High School (Brust Brust and Eschweiler Eschweiler, 
Architects) and the John Audubon High School (Grassold Johnson, Architects). 
some cases the same unit cross section used for several different spans 
order simplify the flasework and detailing. The roofs have been leveled 
with top slabs and the interior spaces used plenum chambers for ventilating. 
The top slabs may either precast poured place. 

Another use folded plate construction long span cantilever hangar 
roofs.(3) The first one built, the Trans World Airlines hangar Kansas 
City, Burns McDonnell, Architects, has cantilever span 150 feet and 
continuous door length over 800 feet. The building 420 feet wide and 
1000 feet long. Preliminary designs have been made for cantilever spans 
250 feet which are perfectly feasible needed. 

second hangar with 150 foot cantilever has been completed for 
World Airlines New York International Airport, and one with 130 foot 
span for Pan American World Airways the same location. Figs. and 
These three have steel bridge cables supporting the outer end. The latter two 
have the folded plate roof curved reduce the positive bending between 
the center building and the cable supports. They all have foot plate 
module poured foot widths which were grouted together after the cables 
were jacked up. order prevent lateral distortion the units which were 
self supporting before being grouted into continuous roof, lateral tie rods 
were cast into the corrugations the top. Fig. shows sixty foot wide 
section the Pan American hangar undergoing load test. This also shows 
the lateral ties and the joint grouted between sections. 

During the load test the TWA hangar, block concrete weighing ten 
tons was accidentally dropped about two feet the outer end the cantilever 
without causing appreciable damage. Both load tests completely justified the 
design assumptions and methods. 

More recently designed are the Mohawk Airlines hangar Utica (Fig. 12) 
with span 120 feet and the National Airlines hangar Miami (Fig. with 
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Fig. Wayne University, McGregor Memorial Community Conference 
Center, Detroit, Michigan (Yamasaki, Leinweber Associates, Architects) 
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Fig. 130-ft. Cantilever Span Pan American World Airways Hangar, 
New York International Airport 


Fig. Load Test Cantilever Roof Span Pan American World Airways 
Hangar, New York International Airport 


111 foot cantilever span. These are built without the suspension cables 
other prestressing. They also have foot module and the folded plate units 
are about feet high the support and three feet high the outer ends. 
all five these hangars, the simple corrugated plate arrangement was 
used with horizontal plates top and bottom between the diagonal plates. The 
plates the two hangars New York International Airport were increased 
one inch thickness over that used Kansas City and recommended the 
designers because ruling the Port New York Authority regarding 
the minimum concrete cover required the New York City building code. 
The American Concrete Institute office building Detroit has cantilever 
roof built tapered folded plates with foot module. The units for each 
side were precast and jointed together after erection welding and grouting. 
Precasting instead pouring place was optional with the contractor. 
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Another interesting application folded plate construction the assembly 
hall now being designed for the University Illinois. (Fig. This will 
400 foot diameter ribbed concrete saucer surmounted dome thin 
tapered concrete folded plates. The saucer will folded plates with top 
slab carrying the seating and the interior corrugation spaces used carry 
heating and ventilating ducts. The rim the saucer will precast con- 
crete ring sufficient strength turn the thrust the dome inward toward 
the central foundation. 

The dome well the saucer will cast segments rotating forms 
reduce the cost. After both the saucer and the dome are concreted, the rim 
will prestressed give them both the proper camber. 

The development the structural design all the preceding examples 
has been the work Ammann Whitney, Consulting Engineers, who also did 
the architectural design where other architect has been credited. 


Fig. 111-ft. Cantilevers 
National Airlines Hangar, Miami International Airport 
(Weed, Russell Johnson, Architects) 


Fig. 400-ft. diameter Dome 
Assembly Hall, University Illinois, Champaign-Urbana, Illinois 
(Harrison Abramovitz, Architects) 
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Methods Analysis 


The following terminology, which consistent with that used Winter 
and Pei(5) will used this paper: Each individual planar element the 
structure called “plate.” The plate “length” the dimension between 
transverse supports, its “width” the transverse dimension between longi- 
tudinal edges and its “thickness” the perpendicular dimension between 
surfaces. The structure “height” the vertical dimension between the upper 
and lower extremes transverse crossection. Bending individual plates 
their planes called “plate action” and bending individual plates out 
their planes called “slab action.” Following list the symbols used: 


Symbol Definition 

Length plate between transverse supports 

Width plate between longitudinal edges 

Plate thickness 

Height folded plate unit 

Distance from support along length plate 

Spacing folded plate unit 

Horizontal projection plate width, 

Angle plate horizontal 

Change angle between adjacent plates 

Load intensity per unit area horizontal projection 

Approximation load intensity (First term Fourier sine 
series) 

Transverse slab reaction due approximate applied load, 

Component slab reaction, plane plate 

Transverse moment plate edge 

Transverse slab reaction due transverse moments, 

Component slab reaction plane plate 

Total transverse plate load (resultant &AS acting 
plate 

Longitudinal shear intensity per unit length along plate edge 

Summation longitudinal shear intensity, 

Moment plate caused transverse load 

Cross sectional area plate 

Section modulus plate for “plate action” 

Moment inertia plate for “plate action” 

Modulus Elasticity 

Deflection plate for “plate action” 
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PLATE CONSTRUCTION 


Symbol Definition 

moments, 

Deflection slab joint caused plate deflection 
Moment inertia unit length plate for “slab action” 
Subscript indicating symbol applicable kth joint plate 
Subscript referring tie between plates 


When shown over symbol indicates value mid span 


The notation taken from Girkmann(4) but has been revised conform 
certain usages common the United States. 


Transversely, the structural action folded plate system similar 
that series continuous frames elastic supports. The sections 
that would isolated series transverse cuts through the structure 
spaced unit distance, may considered elastic frames which are 
supported the intersections the longitudinal plates. The frame reactions 
are supplied the plates acting their own planes. (See Fig. 13) These 
reactions are affected the displacements the plates’ intersections. 


Where the ratio length height small, the “elastic support” 
relatively rigid and the frame action not modified the comparatively 
small deflections. Stress distribution the folded plate can computed 
the equations statics and equation compatible longitudinal stress 
plate intersections. relatively simple method analysis for this condition 
that proposed Winter Pei.(5) relaxation technique, analogous 
moment distribution, used compute the stress variation across the plate 
cross section. Although the method set for simply supported plate sys- 
tems can also applied continuous spans with some modifications. 


For intermediate ratios length height necessary consider 
deflections well compatibility longitudinal stresses. This problem 
can solved either one the following: 


modification the Winter-Pei solution extended include joint 
displacements.(6) The solution involves n-1 simultaneous equations, 
where the number plates. 


equivalent method(7) which involves series simultaneous 
algebraic differential equations including all effects. 


When the ratio length height large, and there series 
repetitions the folded plate unit, the stress variation closely approximated 
treating the transverse cross section that beam, assuming change 
geometry the section. Stresses are then determined standard 
methods beam design. this case plates the end the series should 
investigated for the effects deflection. 

The stress variation can determined using the beam theory for 
length height ratios low when the folded plate unit triangular 
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For Transverse Rigid Frames (b) 
Transverse Unit Frames 
(b) 

Fig. 
increased and their arrangement becomes more complicated, the lower limit 
length height ratio which the ordinary beam theory applicable for 
given configuration cannot completely defined beforehand. However, 
definition the range applicability for specific folded plate system was 
derived Craemer(7) and study his findings indicates that lower limit 
which applicable for the vast majority designs can determined. 
found that this lower limit the length height ratios could expressed 
where constant, dependent upon the folded plate system 
and the remaining terms are indicated Fig. 13. For the ratio 

20, which usual reinforced concrete structures, and for Craemer’s 
configuration, the beam theory applicable single folded plate unit with 
the ratio L/H 9.9. For repetitive folded plate structure the beam theory 


applicable when L/H 7.0. Both these values are much less than the 
ratio L/H which sufficient carry ordinary building loads. Therefore, 
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the use L/H lower limit for structure with symmetrical, repeti- 
tive folded plate units should ensure sufficiently accurate stress patterns 
most building designs. This ratio will lower than necessary for most 
structures. should not applied unsymmetrical units such 
light” structures without further investigation. 

When the ratio L/H the method analysis beam theory out- 
lined below. 


Estimate the unit dead plus live load. 


Isolate transverse strip unit length and determine thickness and 
transverse reinforcing treating the strip continuous beam 
rigidly supported each the slab folds. With this result the esti- 
mated dead load can revised required. 


For one repetition folded plate unit (dimension Fig. 13), compute 
the maximum longitudinal moment. 


For the selected section check the longitudinal concrete flexural stress 
and calculate the required tensile steel, either straight line method 
ultimate strength equations. either event the computation 
similar any *T” beam design. 


Check the web section (diagonal plate) for the effect shear using 
standard shear equations. 


Increase the concrete dimensions required for tensile shear rein- 
forcement. 


Check effect load end plates the discontinuous edge the 
structure determine the adequacy the end plates. This analysis 
will described below. 


For the portion the roof near the end plate analysis similar 
those described previous references(4,6,7) should used. The analysis 
described hereafter adaptation one presented Girkmann, however 
analyses could derived from the other referenced sources. 


Fig. shows end portion folded plate roof Length where the 
roof contains large number folded units. The interior plates this roof 
can designed outlined, but because the lack transverse restraint 
the exterior edge the end plate, there will occur disturbances that ma- 
terially alter stresses and deformations the end plate, becoming smaller 
for successive interior plates. Stress patterns the plates affected the 
edge disturbance are determined isolating the first exterior plates 
the roof. These plates are then treated folded plate roof having zero 
shear force and zero transverse moment the exterior edge plate (free 
end) and zero lateral deflection the interior edge the nth plate. The 
solution becomes more accurate but also more tedious the number plates 
included the analysis increased. has been found that sufficiently ac- 
curate solutions are obtained when the number plates included equal 
1-1/2 units plates Fig. 13, one unit consisting plates). 

The analysis the end plates based continuity stress and de- 
formation along the common edges adjacent plates. each joint, for any 
value (See Fig. 13), the change joint angle caused load, trans- 
verse joint moments and plate deformations and rotations expressed for 
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adjacent plates. Since the joint rigid, the algebraic sum the angle 
changes must zero. For plates possible write n-1 equations 
involving n-1 unknown transverse joint moments. 

the width and thickness the plates and the intensity loading not 
vary with possible greatly simplify the mathematics approximat- 
ing the uniform plate load the first term the Fourier series. The 
approximate applied load 


When equation (1) used, all plate stresses, moments, and deformations are 
expressed sine functions having maximum values mid-span. 

brief, the analysis proceeds follows: 


Assume each plate act simply supported slab spanning between 
adjacent plates and find the loads that would act along the joints (forces 
Fig. 14c) the joints were continuously supported their full 
length. Replace all joint loads their components the planes the 
plates (forces Fig. 14c). Since the plates are much wider than they 
are thick, all longitudinal transfer loads occurs plate action. 


(a) 


(b) 


(c) 


(d) 
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Since the plates are reality continuous across the longitudinal joints 
there will transverse moments the joints that will cause additional 
joint loads (moments Fig. 14d). Compute these additional joint 
loads (AS Fig. 14d) terms unknown transverse moments the 
plate edges. 


Compute longitudinal edge shear forces Fig. 15) for each plate 
required for compatible stresses adjacent plates each joint. 
The transverse loading causing the direct stress are the loads computed 
and 


Compute individual plate deflections caused the transverse loads and 
longitudinal edge shear stresses found previously. These deflections 
will terms load and unknown joint moments. 


Express condition zero angle change (Fig. 16) each joint terms 
applied loads, transverse joint moments and plate deflections. When 
plates are considered there will equations. 


Solve equations from for unknown transverse moments, and then 
solve preceding equations for stresses. 


Fig. cross section through the first six plates the structure 
Fig. taken distance from the near end. Also shown the approxi- 
mate load acting unit length each plate. plate the vertical load 
computed from equation (1). L/2 the load intensity 
Fig. 14c the simply supported transverse slab reaction 
Fig. 14c are the components the slab reactions k-1 the plane the 


adjacent plates. These components are the plate loads. The expressions for 
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These forces can easily obtained graphical resolution force 
the direction the two plates k-1. The positive direction for and 
are measured indicated Fig. 14b. shown Fig. 14d there acts 


along each joint atransverse which has maximum 


value midspan Transverse moments cause changes the resultant 
joint loads the expressions for which are given below. 


(4) 
m -m nm. TX 
The torsional stiffness the first plate neglected and obtained 
treating the first plate cantilever fixed joint The components 
the planes the adjacent plates are computed equations 


AP. 


(5) 


The total transverse load the plane plate obtained summing 
equations and 


the above equations positive when acting away from the joint 
positive when acting toward the edge lower order (i.e., positive 
acting toward edge k-1) 

Fig. shows two successive folded plates isolated along edges k-1 and 
k+1 and distance from the support. and are the moments 


caused transverse loads; and indicate the intensity longi- 


tudinal shears along joint edges; and and are the summation 
these shears from the support section All forces and moments are 
positive shown. joint the direct stress caused all forces and 
moments plate must equal the direct stress caused all forces 
and moments acting plate k+1. This relationship expressed equation 
taken L/2, and become and and the joint shears re- 
quired for compatible direct stress adjacent plates are obtained re- 
arranging terms equation and are given equation 
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PLATE CONSTRUCTION 


Equation and T,, can determined from the assumption 
lateral motion joint Equation written for successive joints con- 
taining n-1 unknown values and unknown values 
The plate deflections, which are obtained twice differentiating the 
expression for are given equation mid span 


solution for unknown joint moments obtained setting all angle 


changes each joint equal zero. Fig. shows the angle changes joint 
caused load, moment, and deformation adjacent plates, the sum 
which must zero. This expressed equation 10. 


— 1 
The angle change caused the load, the angle change caused 
moment and the angle change caused translational deforma- 


tions the plate whole. The expressions for the various terms 
equation are given below. 


(11) 

(13) 


The solution equation (expressed terms eqs. 11, 13) will 
yield values mid span joint moments with which, using equations and 
possible evaluate plate stresses. 
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Equation 10, when written joint n-1 and will require values 
and These values can obtained from analysis the 
typical interior folded plate unit. measure the accuracy the solution 
given comparison the value obtained for from equations and 
that obtained the same joint typical interior unit. 

This analysis requires the solution sets simultaneous equations 
each involving unknowns. The solutions proposed Gaafar(6) when 
applied this problem each require the solution only one set 
taneous equations unknowns. However addition sets “moment 
distributions” must performed evaluate deflection effects and the total 


labor would appear greater than that for the method proposed 
here. 


When the stress the end plate exceeds the permissible value reduc- 
tion stress can obtained, either increasing the size the end plate, 
supplying additional transverse support. The increase the width 
the end plate limited the transverse moment capacity joint Where 
the end plate horizontal, Fig. 14, may possible supply restraint 
against vertical deflection joint the use intermediate columns along 
line The addition these columns would cause reduction thus 
allowing increase bj. 

For the folded plate section shown Fig. the end plate disturbance can 
reduced the introduction intermittent ties between plates and 


single tie were added mid span the change deflection plates and 
given below. 


(14) 

Ast (15) 


The force the tie obtained from equation 16. 

the above equations the force the tie and and are 
the crossectional area, length and modulus elasticity the tie respectively. 
With the addition the tie approximate solution for plate stress ob- 
tained substituting Ait and for the original values 


deflection equation 13. revised solution for joint moments and plate 
stresses can then obtained. 


Standards and Practice 


Except for few provisions particularly pertinent folded plate construc- 
tion, practice connection with aggregates, material and product, tests, 

concrete quality, allowable unit stresses, the mixing and placing concrete, 

forms and details construction, etc., should follow the normal require- 


ments codes such the Building Code Requirements for Reinforced Con- 
crete (ACI 318-56). 
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Folded plate design will vary from what specified the typical code 
regard methods analysis, precautions against buckling instability, cam- 
ber and deflection controls, permissible combined flexural, shear and axial 
unit stresses, corrosion control and fire protection requirements. 

Folded plates, like shells, utilize thin sections with the concrete distributed 
over large surfaces. The load carrying capabilities these sections, be- 
cause thinness, depend the local well the general stability the 
section. Readily available methods analysis for buckling stability plates 
with varying types edge support may used. Somewhat higher factors 
safety are recommended than normally required metal construction. 
most concrete folded plates factor safety can provided 
small variations the thickness and geometry the folds without affecting 
the economy the structure. Structures designed this criteria have be- 
haved well under severe overload and test conditions. Designs based upon 
ultimate strengths for safety under overload should also consider the local 
effects stress transfer ascertain that the assumed configuration will 
either maintained ultimate loads or, distortion occurs, that the re- 
duced capacity adequate provide the necessary factor safety. Changes 
stress caused geometrical displacements should included the 
design. 

Camber and deflection controls follow those prevalent for beam and slab 
design. Overall deflection can generally determined the basis the 
gross crossectional area for computing the moment inertia. For longtime 
loads the effective modulus may approximated using 750,000 psi for 
stone concrete and 500,000 psi for lightweight aggregate concrete. Where 
sag becomes problem, small amount compression steel will very 
effective controlling deflections. Other local deflection additives may 
approximated consideration the effects transverse distortions 
described the methods analysis. case the plate width narrow 
compared the span and the break the surface the connecting edges 
small, strains the plane the plate may cause large local displacements, 
loss edge support and consequent buckling the member. Rules account- 
ing for the latter motions are difficult formulate because the infinite 
combinations possible different configurations, but for angles over 25° 
plates 1/20 more the span, the effect may considered small. 

Care must exercised not introduce curves into tension chords 
during the cambering process the introduced radial force components will 
tend force the edges into straight lines. Such cambers are useless for 
correcting overall elastic and plastic deflections. 

Unit stresses can governed usual code requirements, however, care 
must exercised relate maximum stress buckling stress adequate 
factors safety. Except where governed buckling, direct design stresses 
the flange and shearing stresses (with web reinforcement) 


200 psi have been used numerous times without signs distress. 

The usual practices relating corrosion protection are applicable 
folded plate construction. While the cracking the deep corrugated sections 
may less than slender conventional beams, the light reinforcement 
somewhat more difficult hold fine tolerances, particularly the sloped 
sections and there appears justification for reducing cover required the 
ACI code buildings comparable construction. Carefully controlled pre- 
cast folded plates with accurately placed reinforcing steel and rich concrete 
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steel coverage should governed the presently established ACI require- 
ments for Thin Skin Precast Concrete. 

Numerous inadvertent tests have proven that folded plate structures with 
the usual ACI covers behave well under severe fires. Intensities high enough 
spall and strip the cover over the bars have not caused collapse, disturb- 
ing deflections or, demonstrated subsequent load test, reduction the 
safe capacity the member. Assuming the use suitable aggregate, cover 
sufficient satisfy corrosion protection requirements should sufficient 
for fire protection. Where fire causes spalling, subsequent protection 
reinforcing steel can provided use gunite patching other suitable 
remedial measures. 

points restraint such stiffener ribs, frames and end beams special 
consideration should given the plate design. Ordinarily the folded plates 
are able accommodate volumetric changes slight sidesways and small 
flexural adjustments. However these thin plates are subject rather quick 
and intense volume changes and the sections invariably will crack where 
these motions are restrained thicker slower acting frames and rigid 
beams. these restraint points, the slab must reinforced the juncture 
resist and distribute cracking. 

Both lightweight and stone concrete have been used extensively for folded 
plate roofs. Both provide the strengths required for any normal design. 
While the modulus the lightweight concretes may only one half two 
thirds that stone concrete, this usually does not cause deflections which 
cannot provided for the design. 


Construction Procedures 


Formwork and Falsework Forming details for folded plate construction 
have been varied the applications that can made folded plates. 
majority the simple span falseworks have consisted light wood 


tubular working platform which set self-contained surface units shaped 


the fold configuration. Space may provided the platform that the 
upper sections can dropped integral unit means jack assembly 
located ground level. one installation small construction holes were left 


the folded shell. When the concrete had obtained sufficient strength the top 


sections were supported from the shell through these holes while the lower 
platform was moved forward, after which the top sections were lowered 
ground level means cranes and subsequently rolled forward that 
level. This procedure added the number operations but obviated more 
complicated raising and lowering operations either the platform top, 


ground levels. 


large, deep folded plates the form sections are more intricate and 


heavier and usually necessary provide mobile falsework accommo- 


dating the greater weights and the deeper sections.(8) Various schemes have 
worked satisfactorily and doubtful whether even the contractors have 
sufficient cost data evaluate any appreciable savings one scheme over 


the others. 


Because the desirability obtaining maximum re-use the falsework 


without excessive delays progress the work customary permit 
decentering the formwork within days from the day the pour. 


Decentering controlled stripping cylinders which are required show 
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minimum specified strengths prior the decentering. These strengths ordi- 
narily are established assure minimum modulus and strength least 
2500 psi the time stripping. 

Decentering the falsework ordinarily does not invoke particular diffi- 
culties but requires solely that the steps the operation are recognized. The 
falsework must designed that strains caused any part the pour will 
not damage concrete already reaching the set stage. The decentering can 
usually done using visual reference line level points the concrete 
structure will usually much stiffer than the supporting falsework. this 
time the concern limited obtaining gradual easing individual post 
loads such that excessive load concentrations are not induced individual 
framing and supports which could overload the falsework locally, buckle in- 
dividual posts, provide reactions the plate structure that were not con- 
sidered the design. 

Joint and Edge Details Construction joints the edges the plate pours 
usually are not critical because the flexibility the edges usually permits 
closure succeeding pours shoring clamping adjacent edges into one 
plane. One common method making the joints decenter several sec- 
tions and then connect the adjacent sections means narrow cast joint. 
some cases, where falsework area available, the edge the previous 
pour shored the theoretical position and the new concrete cast against 
the previous section. The latter method will result load transfer between 
sections but the folded plate unit usually too flexible bending and torsion 
reduce the ultimate load capacity the sections. 

The edge configuration may handled many ways depending the span 
the folded sections and the desired architectural treatment. many types 
folded plate structures vertical horizontal edge plate will effectively 
terminate the sides the building. For large spans the width horizontal 
plate needed restrain the edge the adjacent plate will require thick sec- 
tions and heavy reinforcement carry its own weight cantilever. this 
case often advisable cut down the required width tieing the exterior 
plate back interior plates and thus reducing the “plate action” bending. The 
temporary edges construction joints can treated the same manner. 
Columns may placed under the edges for vertical support properly con- 
sidered the 

Reinforcing Reinforcing steel likely light and closely spaced 
distribute such cracking may occur the concrete. Because this light- 
ness and the difficulties added the sloped plates, the steel place should 
protected walkways and working platforms. The steel patterns are 
likely appear complicated the drawings but only because the steel 


usually more elaborately detailed than for more conventional slab construction. 


Folded plate structures are usually long with simple steel patterns. This 
results arrangements ideally adaptable ground fabrication much the 


heavier steel into mats. These mats can placed large sections means 


crane operations. such cases careful attention should given the lo- 


cation each layer bars that the fabricated systems will mesh properly 


all points without lacing the bars. 

Casting Techniques Casting techniques for folded plate construction are 
not particularly difficult different from that for conventional construction. 
The overall section usually effective that the slabs are insensitive 
local flaws weaknesses except possibly the areas maximum shear 
transfer the support points. the transverse direction between folds, the 
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strength requirements are similar those normal slab construction except 
that the span usually short that the flexural stresses are low. Accord- 
ingly the folded plates can cast employing the same concrete handling tech- 
niques utilized for normal concrete joist, waffle and flat slab construction, 
spite the thinness construction and some difficulties added the sloped 
surfaces. 

the slope the surfaces approaches 45° care must exercised 
achieve the desirable compactness and smoothness. the experience the 
writers most sloping plates 45° have been poured without the use top 
forms. accomplish this, the bottom plates are usually placed ahead the 
rest the pour. This followed feeding the concrete down the slope 
from the top and working the concrete around horizontal bars the slope. 
Concrete slumps maintained between 1-1/2 2-1/2 inches are stiff enough 
for slopes 45°. Higher slumps than these will result flow the con- 
crete down the slope while lower slumps are unnecessary and too unworkable 
allow reasonable production speed. Bars closely spaced, which desirable 
for stress distribution and reduced plate thickness, will also support the wet 
concrete and help prevent from slumping down the slope. Spading 
preferable excessive vibration which may initiate sag the finished con- 
crete immediately behind the pour operation. 

other projects, methods pouring sloping surfaces have included the 
use modified top forms. Examples these range from light ply-wood 
panels attached the lower surface large size crane-handled sectional 
frames covered open mesh surface. These double forming methods 
allow vibration compaction, the open mesh being particularly effective 
restraining the concrete while permitting escape entrapped air and water. 
Methods such these utilizing-the lower form surface for alignment and 
support are inexpensive, running from $0.05 $0.12 psf. 

more general method pouring utilizes rack ladder section with 
longitudinal dividers set above the concrete surfaces (Fig. 17). The section 
progresses with the pour and forms ladder-platform for the workmen which 
also restrains the concrete from slip flow during the vibration. 

new job, the techniques concrete control, slump regulation and 
vibration are usually mastered the first three four hours the first 
pour, after which the pouring operation proceeds smoothly other 
construction. 


Protective Cover and Insulation 


wide variety covers and insulations have been used folded plate 
structures with varying degrees success. Among the coverings are in- 
cluded reinforced bitumens, cold bitumens, bituminous roofing felts, asphalt 
and reinforced asphalt, vinyls, acryllics and acryllic plasters. Unfortunately 
the experience record too short most these roofing types offer 
more than comment and such recommendations are believed sound 
the basis limited experience. 

The folded reinforced concrete sections are inherently tight waterproof 
type construction. Usually the supporting structure these roofs offers 
little restraint the direction the long span and the folded slab flexible 
the transverse direction they are not subject the usual shrinkage and 
temperature cracking inherent restrained slab construction. When 
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Fig. Placing Concrete 


prestressed the main span direction and reinforced distribute local 
cracking stiffener restraint points, the shells are watertight and have been 
tested six foot water heads without leakage. With non-prestressed 
reinforcing, minor cracking will occur the bottom plates simple spans 
and also along the top plates points restraint cantilevered areas. 
Again these cracks are limited the reinforcing well distributed and 
worked normal unit stresses. However, spite the good experience 
with bare plate roofs, the practice leaving the concrete devoid insulation 
and covering invites larger motions and water penetration and highly ques- 
tionable where the class construction warrants cover. 

Normal multiple felt roofing laid over insulation and anchored wood 
strips will assure protection and can furnished with 10, and year 
roofing bonds. The cost will range between $0.30 $0.50 psf. 

Vinyls are available variety rich colors suitable for exciting archi- 
tectural effects. Experience has indicated that these vinyls will behave satis- 
factorily except over wide working cracks that open after the coating 
applied. Such cracks, incidently, can covered vinyl reinforced with 
glass fabric. However, vinyls form membrane sealing the surface the 
concrete with the possible entrapment moisture and consequent blister 
problems. Adequate thicknesses (30 mils) must used and the cost 
will likely range from $0.50 $0.75 psf. 

The experience with acrylics has been both good and bad depending the 
circumstances. Placed under specified conditions clean surface, the 
bond has been good observed over several years and this material can 
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supply acceptable cheap roof. Placed over wet concrete areas where 
water ponds, dusty surface applied during low temperatures, this cover- 
ing has been known strip badly over large areas. The colors available are 
good. Solvent type and water emulsion types are available with the advantage 
appearing with water emulsion because the probable moisture content 
the concrete. The acrylic coatings will not bridge cracks formed after ap- 
plication the surfacing. The cost very low, being only about third 
felt ply roofing. Acrylics reinforced with glass fabric may good cheap 
roofing solution but experience record available such installation. 

good solution where insulation not required above the slab surface 
appears the use asphaltic coating with powdered aluminum surface 
and reinforced with glass fabric. The reinforcing flexible enough bridge 
small cracks. Aluminum green, blue and limited red colors are available 
with all apparently reducing aluminum The cost somewhat less 
($0.20 $0.35 psf) than that ply roofs. Ten year guarantees lower 
cost and year bond higher cost will furnished certain manu- 
facturers. 

Insulation materials have been applied both the inside and outside 
folded plate slabs. Exterior insulation desirable that reduces 
volumetric strains the concrete and keeps the dew point temperature out- 
side the concrete surface. Such insulation usually set between wood 
nailers and covered with felt roofing. 

Insulation the undersurface the concrete has been used many 
installations with reported good results. many cases this inside insulation 
has served the double purpose acting the form surface dual purpose 
insulation and acoustical material. Inner insulations permit the use 
roofing membranes offering wider selection colors. 


Economics 


Basic data typical material quantities are presented Figs. and 19. 
This information based particular folded plate units, section heights, 
control joints and type beam framing. The expansion this information 
estimated building costs, shown Fig. 20, furnishes interesting com- 
parison with joist construction. These serve illustrate the possibilities 
folded plate construction. Fig. also reflects contractors’ bid prices for 
specific designs built widely separate locations utilizing varied construction 
methods. 

The average quantity concrete required per square foot floor area for 
folded plate structures will obviously depend the configuration the folded 
plate utilized the design. Average thicknesses based some prelimi- 
nary designs and checked against number final designs are shown Fig. 
18. These designs are largely for plates folded into corrugated form simi- 
lar shape that the Central Warehouse roof. (Fig. 10) For simple 
spans the ratio span height ranges from 25. For cantilever spans 
the ratio span height approximately 10. 

For short spans the thickness the folded plate controlled the cri- 
teria established for protection reinforcement. For the designs used 
cost comparisons minimum cover 3/4 in. transverse steel and 1-1/2 
inches main steel was selected. Minimum steel and these covers governed 
the plate thickness spans less than about ft. where flexural and shear 
stress requirements begin govern the design. 
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AVERAGE THICKNESS CONCRETE INCHES 


SPAN FEET 


Fig. 


was assumed that the folded plates terminate continuous beams 
each end each span. Inasmuch these beams vary only slightly total 
dimension with increases the span, the average thickness (of slab plus 
edgebeam) actually decreases with increase span optimum point. 
This variation reflected the plot the average thickness curves. 

the span increases beyond the point minimum materials for the 
given sectional height, lesser quantity materials could obtained than 
shown the curves decreasing the span-to-height ratio. However the 
entire configuration the plates used any given project controlled 
many other factors such esthetics, building volume, roofing and pour- 
ing costs, etc., that refinements were considered beyond the scope this 
paper. 

The reinforcing steel requirements are relatively low for folded plate 
construction, Fig. 19. This results from the unusually large lever arm fur- 
nished the folds the direction the main span compared flat slab 
concrete joist construction, and the usually very short transverse span 
between foldedges. Furthermore the bulk the concrete dispersed into 


areas and geometrical positions where highly efficient carrying the 
developed stresses. 
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The cost concrete placement similar that for flat slab and joist 
construction with the exception some additional time required for the 
handling and compaction the concrete the inclined surfaces. ex- 
ample, the first and always most difficult pour ft. high, ft. span, 

pour section covering 3200 square feet roof area folded plate and edge- 
beams used c.y. concrete. crew consisting one foreman, ten 
laborers, and two cement finishers placed and finished all work the roof 

7-1/2 hours. crane with operator and two mixers and four man crew 
mixed and delivered the material the roof. 

Cost estimates for folded plates corrugated configuration are shown 
and compared with similar estimates for concrete joist construction Fig. 20. 
The unit prices used the cost buildup are follows: 


Joist Folded Plate 
Concrete $22/Yd. $37/Yd. 
Reinforcing $0.16/# 


Formwork $0.60/SF $0.85/SF 
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Fig. 


These estimates show that generally the conventional framing practical 
spans 40-ft. Beyond this the cost increases rapidly and diffi- 
culties are introduced due weight and deflections which limit further in- 
crease span. The folded plate spans however can increase ft. 
more cost only $0.40 $0.50 per square foot over the cost the 
short spans utilizing conventional framing. This true for both lightly 
loaded roofs and heavily loaded floors. The flat surfaces and geometrical 
form many cases permit utilization the exposed concrete surface for 
architectural treatment. This results the elimination suspended ceil- 
ings and savings thus effected may more than the increase the long span 
folded plate construction compared with short conventional spans. 

Another saving achieved certain installations has resulted from the 
utilization closed cells ducts distribute conditioned air. The savings 
shown contractors’ prices for alternate schemes with and without metal 
ductwork one particular project showed savings great enough pay for 
substantial part the total cost the roof. The results experienced this 
building over several years operation utilizing the structural system the 
duct lines were exceptionally good because the area and volume the ducts 
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were larger than could afforded normal mechanical installations. This 
resulted low pressures, low velocities, and draft-free rooms with almost 
uniform temperature throughout. Estimates have indicated the economy 
this framing special applications. This illustrated one installation 
involving ft. spans carrying 250 psf live load (Fig. 11). Here numerous 
preliminary designs structural steel, concrete beam and slab construction, 
and composite structural steel and concrete emphasized the difficulty 
building structure which would allow passage the large ductwork systems 
needed for air conditioning within ft. total height. ft. high closed 
folded plate ft. module was designed and constructed without difficulty. 

Where folded plates are used roof construction, care must exercised 
shape the roofs that drainage can accomplished without multitudinous 
headers relieve local trapped areas. this not done the added plumbing 
costs will tend price the overall structure out competition. 
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LATERALLY DEFLECTED COLUMNS 


John Sherman,! ASCE 


SYNOPSIS 


Column formulae current are based the assumption 
that the unbent axis the column remains its original position for all con- 
ditions loading. the object this paper consider the case 
column which the unbent axis rotated one end due the deflection 
the other end, condition which occurs portals and bents when subjected 
lateral loading, and truss members due the deflection the truss. 
For the purpose this paper column which deflects described above will 
called laterally deflected column. 

The simplest case laterally deflected column has both ends fixed 
prevent roration. When the column loaded with axial compressive force 
and horizontal shearing force applied the top the column, the 
deflection the column can considered made two parts, the de- 
flection due and the deflection due the displacement the force 
The deflection due produce additional bending which increases the de- 
flection, this phenomena continuing until the deflection reaches final value. 
will shown how the value this final deflection can determined from 
the differential equation the elastic curve the column. 

The moments the ends the column, which are equal for the fixed end 
column, can closely approximated follows. the deflection due 
ignored and the deflection due then the bending moments the 


ends the column become 


which the bending moment and the length the column. (1) 


forms the basis the approximate method which will discussed this 
paper. 


Note: Discussion open until March 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2221 
part the copyrighted Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. October, 1959. 


Engr., American Bridge Div. Steel Corp., Roanoke, Va. 
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NOTATION 


The letter symbols adopted for use this paper are defined where they 
first appear, the text illustrations, and are arranged alphabetically for 
convenience reference the appendix. 


INTRODUCTION 


Column action has long been matter great interest and the present 
formulae used design closely approximate the results obtained from mathe- 
matical investigations and loading tests. For the case where the axial com- 
pressive force eccentrically applied the ends the column the secant 
formula applicable. 

Considering the axial compressive force load, the secant formula 
covers the following conditions application the load 


The eccentricity application the load the same each end the 
column. 


(a) The eccentricities lie the same side the column. 
(b) The eccentricities lie opposite sides the column. 


The eccentricity application the load larger one end the 
column than the other end. 


(a) The eccentricities lie the same side the column. This the 
general case the secant formula which can stated follows; 
there magnitude the load which will denoted which 
will produce moment point the column, near the end where 
the large eccentricity applied, equal the moment that end 
the column. Therefore, any load larger than will produce 
larger moment interior point the column than the moment 
the end where the large eccentricity applied. 

(b) The eccentricities lie opposite sides the column. this case 
the maximum moment occurs the end which has the largest 
eccentricity. 


The column subjected lateral loads addition the load 
this case the eccentricities are taken equal and the same side 
the column. 


all the above cases the secant formula the unbent axis the column 
assumed remain position, therefore, the secant unless modi- 
fied not applicable the columns the laterally loaded portals shown 
Fig. 

The column bases Fig. are fixed, and the portal deflects the columns 
rotate about the bottom panel point the portal relieving the moment this 
point and increasing the moment the base the columns. 

the location the point contraflexure and the deflection this point 
are calculated from the moments due the shearing force then approxi- 
mation the moment the base the column 


HmL+ (2) 


lat 
eac 


flex 
Fig 
q 
the 
des 
fec 
cor 
slo 
pre 
des 


ere 


ase 


odi- 


lumns 
this 


point 


(2) 


COLUMNS 


Fig. Showing Typical Portals. 


which the distance from the base the column the point contra- 
flexure. 
The portal shown Fig. 1(b) has shorter panels than the portal shown 
Fig. 1(a) and has smaller moments the bases the columns. Therefore, 
using stiff portals the effect lateral buckling can controlled design. 
This paper will restricted the analysis columns such shown 
Fig. for which the following assumptions will made; 


(a) The unit stresses the portal bracing small and the rotation the 
portal can neglected. Therefore, the panel points the portal will 
assumed remain vertical line during deflection. 

(b) Since the columns are much more rigid than the portal bracing will 
assumed that the column hinged the top the portal and con- 
tinuous over the portal struts. 

(c) From the above assumptions follows that for all loads applied within 
the depth the portal the shearing force and the moments and 
are the same for each column. 


The above describes the arrangement and loading for which equations for 
the safe working stress will calculated, the approximate method 


described, and the more precise method which takes into account the ef- 


fect the deflection due the displacement the line action the axial 


compressive force 


The case column fixed both ends will first considered: then the 


slope the column the bottom the portal and the effect the axial com- 


pressive force this slope will calculated; and finally the two cases 
described above will combined cover the general case shown Fig. 


The Column Below the Portal 


Fig. 2(a) shown simply supported beam having length with 
lateral load acting midspan and axial compressive force applied 
each end. 
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Simple Beam Laterally Deflected Column 
Fig. Showing Elastic Curves. 
Fig. 2(b) shown laterally deflected column fixed both ends with 
lateral shearing force and axial compressive force applied each end. 
The beam Fig. 2(b) the semi inverted form the beam Fig. 2(a) 
and has the following similarity; 
flexu 
moment distance from the ends the laterally deflected 
column numerically equal the moment distance from the 
midspan the simple beam, therefore, the moments acting the ends 
the laterally deflected column are numerically equal the moment 
midspan the simple beam. for 
(b) The curvature any point the laterally deflected column numeri- 
cally equal the curvature the same point the simple beam. 
(c) From the above follows that the deflection the point contraflexure 
the laterally deflected column equal the deflection midspan 
the simple beam. 
seri 
The solution for the moment midspan the simple beam known, (2) 
from which follows that the moments the ends the laterally deflected 
column are 
which tern 
tern 
vel (3a) sma 
Eq. (3a) the radius gyration the plane bending, the aver- 
age compressive stress due and the modulus Elasticity 
and compression. 
Eq. (3a) can written 
Safe 
from which 
Substituting for Eq. (3c) the value the critical load 


COLUMNS 


This the smallest load that can keep the column the deflected shape 
shown Fig. 2(b). The load the Euler load for the column, the de- 
flection due this load apparently indeterminate. This because the ap- 


proximate expression for curvature has been used deriving Eq. (3). 
The value the average compressive stress obtained from Eq. (3a) 
(5) 
which when taken becomes the buckling stress fp, where 


a 


When the deflection due ignored the deflection the point contra- 


flexure the laterally deflected column 


(3c) 


(4) 


7 
7 


(7) 


Substituting Eqs. (3c) and Eq. (2) gives the following approximation 
for the moments the ends the column 


The difference between Eq. (2) and (8) shown the following infinite 
series 


which seen that the first two terms appear Eq. (8), the remaining 
terms representing the effect When equal taking the first two 


terms the series the error six and one half per cent. Therefore for 
smaller values Eqs. (2) and (8) are good agreement. 


Equating the yield point stress the total stress the column 


(10) 


which the allowable average compressive stress based using 
safety factor the distance from the neutral axis the column the 


4 
7 


extreme fibre, and 


(12) 


Solving Eq. (10) for 


4 
4 
x 
= a 
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lexure 
(2) 
(3) 
aver- 
ension 


Values are given the appendix. 


The corresponding equation for based the approximate method 


has been stated that the secant formula not applicable lateral de- 
flected columns unless modified. The following form the secant formula 
noted Formula AASHO specifications. 


(14) 


the above formula, the bending moment due loads applied trans- 
verse the axis the column, the deflection due the ec- 
centricity the axial compressive force which taken act the same 
side the column each end, corresponds fy, sec corresponds 
and the factor .25 taken represent the inherent crookedness the 
column; the remaining terms correspond the notation used this paper. 

This formula approximate the same extent Eq. (12) except for the 
secant term which additional approximation since correct only for 
the case continuous uniform load. 

Taking the transverse load the column can semi-inverted and can 
taken apply the laterally deflected column making the following 
substitutions; and using Eqs. (3c) and (8), from which the 
following equation finally obtained. 


than one. 


The factor safety taken 1.76 AASHO specifications, using this 


value for the Euler stress becomes equal when 


33,000 


Therefore, the yield point stress 33,000 p.s.i. can used for Eqs. 
(11),(12) and (15) when applied primary members. 
Eqs. (11) and (15) must solved trial and error, Eq. (15) giving results 
which compare favorably with the results obtained Eq. (11) but more 
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troublesome solve. Eq. (13) quadratic respect and can used 
obtain trial value use the solution Eq. (11) (15). 


The Column Within the Portal 


the following general equation, which takes into account the effect 
the axial compressive force will developed for slope the elastic curve 
(13) the column the bottom the portal. 

The situation the column within two panel portal shown Fig. 
de- axial compressive force applied each end, and moment ap- 
ula plied support which represents the bottom the portal. can due 

any loads the column below the portal the portal panel points. 

The slope 69, support can determined when the relation between 


Substituting Eq. (17) the differential equation the elastic curve the 


following differential equation the second order obtained; 
the 
The complete solution Eq. (17a) 
the 


The constants integration, and are determined from the condition 
that zero each end spanaL. is, therefore, concluded that 


and 
B=- 
Ptankal 


this Substituting the above values and (17b) gives the following 
equation for 


(17e) 


Eqs. (18) 
results 
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Fig. Showing Elastic Curve Column Within Portal. 
The first derivative Eq. (18) 

the slopes supports and 


(20) 
Va 


which 


(20b) 


kal 


Using the notations pug and Eqs. (19a) and (19b) take the following 
abbreviated forms, 


(21) 


The factors pug and are found tabulated existing argu: 


ment 2u, and represent the effect the axial compressive force the 
moments and and, therefore, the effect the slope the elastic curve 


whi 


loa 
tha 
(24 


609 
wh: 
the 
Introducing the notations 
anc 
q 


the column the bottom the portal. For any given value the load 
increased the values and increase very rapidly. For ex- 
ample when increased from .50 1.00 and increase 35% and 
60% respectively. 


From the condition that support assumed hinged the slope the 

elastic curve the left end span found 

where 
(22a) 


Substituting Eq. (22) Eq. (21b), the following relation between and 


obtained 


which when substituted Eq. (21) gives the following equation for the slope 
the elastic curve the column the bottom the portal. 


(23) 


(19a) 
(19b) 
Denoting the bracketed term the right side Eq. (24) where sub- 
script denotes the number panels the portal, Eq. (24) becc nes 
(20) which for one panel portal 
apy 
(20a) 
and for two panel portal 
(20b) 
When the effect ignored, then 
(24e) 
(21) 
For portal having any number equal panels and the effect the axial 
(21a) not included, the moment the second panel point can not larger 
(24e) .875a. For infinite number panels the value .867. From 
the which seen additional stiffness gained using more than two panels 
curv 
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portal. However, diminished the column tends fixed the 
bottom the portal, providing the portal has sufficient depth resist 
rotation. 

From the above seen that proper design can keep the effect lateral 
deflection the column minimum. 


Combined Action the Column and Portal 


Fig. shows the situation the column below the portal; the column 
fixed the base and rotates through angle the bottom the portal, 
which deflects distance reference the tangent the elastic curve 
the column base. 

axial compressive force applied each end the column and 
lateral load applied distance from the base the column, which 
are resisted lateral shearing force applied the bottom the portal 
and the moments and applied shown. 

The equations for the bending moment the left and right the load and 
the bottom the portal are 


(25a) 


Substituting Eqs. (25) and (25a) the differential equation the elastic 
curve then for less than 


and for larger than 


When Eq. (26) and its first derivative are zero, from which con- 
cluded that 


(26b) 


(Q-H) 


the point application the load the two portions the elastic curve 


given Eqs. (26) and (26a) have the same deflection and common 
then 
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ral 

Fig. Showing Elastic Curve Column below Portal. 


(25) from which 


(25a) (27b) 

(25b) 

tic 


Substituting the above values and Eq. (26a), then when 


(28) 
(26a) 
Substituting the above value Eq. (25b) the following relation between 


When the first derivative Eq. (26a) equal which given 
Eq. (24a), then 
(26c) 
Substituting the right hand side Eq. (29) Eq. (30) gives the following 


for 
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(31) 


Substituting wdcL for QcL Eq. (31) and integrating between the limits 
zero and gives the following equation which added the first term the 
right hand side Eq. (31) will give the value due continuous uni- 
formly distributed load having intensity extending from the base the 
column the bottom the portal 


When the only lateral load applied the column Eq. (31) becomes 


the denominator the second quadrant, and negative when 


larger than and the value negative when infinite. table 
containing values for the second quadrant the appendix. 
When zero the column fixed the bottom the portal and Eq. 


(31b) Eq. (3). 


Values are shown Fig. argument for several stated 
values 


The heavy vertical lines Fig. are the asymptotic values 


which becomes infinite for equal 1.00, .80, .60, .40 and .20 
respectively. 


The buckling stress occurs the elastic range for values larger than 
given the following equation, 
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Values 
(32) 
r 
Tabulated below are the values beyond which buckling occurs the 
elastic range. 
Eq. noted from Fig. that increased the column approaches 
hinged condition the bottom the portal, from which concluded that 
ated the distance the point contraflexure the column then 
(33) 


decreased approaches unity, therefore, when small the 


effect the axial compressive force ignored, then can approxi- 


mated where 


(33a) 
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Using (31c) the formula for the safe working stress when portal action 
taken into account becomes 


the approximate Eq. (12), and the secant formula Eq. (15) the approxi- 
mation for not linear and the correction for portal action can not 


applied Eq. (34). apply the axial compressive force the point 
contraflexure Eq. (8) must modified follows 


(35) 
and with this modification the approximate and secant formula become 
(36) 


(37) 


modifying the secant formula cover the case lateral deflection the 
factors and .25 were omitted, for the reasons that follow. 

Since known calculated eccentricity the stresses calculated from 
the moment are secondary. The stresses calculated from the moments 
portal are also secondary, and since these moments are calculated there 
need for the factor the equations for the working stress. 

The following statements can found Steel Column Research Final 
Report the Special Committee. (5) 


(a) The factor .25 assumed contain allowance 0.15 for secondary 
stress and 0.10 for crookedness. 

(b) For members bent reverse curvature the effect crookedness 
small and need scarcely considered. 


AASHO Specifications the secondary stress represented and .25 
taken represent the effect inherent crookedness regardless the 
nature the curvature. 


the following comparison will made the effect inherent crooked- 


ness pin end column and laterally deflected column. 

Assuming the shape pin end column when the loading removed 
sine curve with mid-ordinate has been that when load 
applied the ends the column the equation the elastic curve 


and 
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(38) 


which 


(38a) 
For laterally deflected column with fixed ends when the loading re- 


moved the mid ordinate would and the shape the crooked column could 


described the following equation 


which satisfies the condition that and are zero when and 


Placing Eq. (39) the differential equation the elastic curve for later- 
ally deflected column with fixed ends, and integrating, then 


which 


(40a) 
When and are zero, then 


(40b) 


Substituting Eqs. (40b) and (40c) Eq. (40), and adding Eq. (39) gives the 
following equation for 


and 


Differentiating Eq. (41) and letting found that the moment 
tanu 


(42) 
From Eq. (38) the maximum moment mid-span the pin end column 
(43) 


T8 
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comparison the effect inherent crookedness for the pin end column 
and the laterally deflected column given below. 


Equation 1.22 1.33 2.22 


From the above seen that the effect inherent crookedness need not 
considered for laterally deflected columns. 

the equations this paper length has been used which has not been 
qualified. The recommendations for reduced length AASHO specifications 
not apply. laterally deflected column has reverse curvature and Steel 
Column Research Final Report The Special Committee stated that 


little consideration given the length columns with reverse 


ture the maximum moment occurs the ends. 

However, portal column the column becomes reinforced enters 
the gusset and would conservative measure the length the column 
point midway between the top the floor beam and the upper edge the 
gussets. 

Denoting the three formulae for the working stress laterally loaded 


column the Tangent Formula, the modified Secant Formula and the Approxi- 


mate Formula, comparison the stresses calculated these three formu- 


lae will given for silicon steel column. will taken and 120, 
1.00, .50 and zero, and .02, .06 and .10. 


For and 120 buckling occurs the elastic range for 1.00, 


and 1.00, .50 and zero, respectively, and the unit stress use the de- 


nominator the equations for the working stress obtained from Eq. (32). 
the tabulation below the safe working stresses are given kips per 
square inch. 


Tangent Formula Modified Secant Formula Approximate Formula 
4.2 3.9 8.0 4.7 
3.9 2.3 4.2 2.5 
2.8 1.8 2.8 1.9 
4.5 9.0 4.4 10.3 5.6 
5.0 2.6 5.7 3.1 
3.7 3.7 2.0 3.8 2.3 
10.2 5.0 10.0 4.8 11.8 6.9 
5.9 2.9 5.8 2.9 6.0 3.4 
4.2 2.3 4.3 2.8 
CONC LUSIONS 


has been shown that the modified secant formula gives results which 
compare favorably with the results obtained the tangent formula. 
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deflection not contained A.R.E.A. specifications, however portals rail- 


flection. There column formula specifications which can 


belong the class laterally deflected columns. The formulae this paper 
can applied such truss members after determined from the end 
moments. 


stresses this paper, seen that the modified secant formu- 


based applying the axial compressive force the point contraflexure 
from Therefore, the consideration buckling limited de- 


use design. fact this the only practical method use when the 
moment inertia the column variable. Furthermore, since the bending 
stresses are secondary nature refined calculations should not required. 


occurring erection cantilever bridges, the precise method should ap- 


Therefore, AASHO Specifications contains formula which properly handled 
applicable laterally deflected columns. 
formula the type which can modified cover the case lateral 


way bridges are invariably stiff that lateral deflection need scarcely 
considered. 


building construction long columns with shallow portals often occur the 
minimum depth which can determined calculating the lateral de- 


modified cover the case lateral deflections. 
Truss members having reverse curvature have their ends displaced and 


For the proportions portals, and the ratio 


used computing the 


sufficiently precise. This formula and the approximate formula are 


termining the stress replace the denominator these formulae. 

The approximate formula not precise the modified secant formula, 
particularly when the stresses are small. Since the conditions the ends 
the column are not definitely known this formula should satisfactory for 


When the column very slender and the portal deep, condition frequently 


plied obtain the safe working stresses. 

The precise method, from which the tangent formula was derived, makes 
the usual approximation for curvature. Therefore the angle approaches 
its critical value the method becomes less precise and should not used. 


APPENDIX 
Notation 


Modulus elasticity tension and compression. 
Horizontal shearing force. 
Moment inertia the plane bending. 


Precise magnification factor for portal action. 
Approximate magnification factor for portal action. 
Length. 


Bending moment. 


Axial compressive force 


Concentrated load. 
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factor the slope the elastic curve the column 
the bottom the portal. 


aL, bLandcL Lengths. 
Distance from the neutral axis the extreme fibre com- 
pression. 
Yield point stress. 
Distance from base column the point contraflexure. 
_kL 
1.0 
Deflection. 
Slope elastic curve. 
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1.000 
1.003 
1.014 
1.031 
1.057 
1.093 
1.140 
1.204 
1.287 
1.400 
1.557 
1.786 
2.144 
2.771 
4.141 
9.400 
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TABLE 
tanu First Quadrant 


1.000 1.001 1.001 1.001 1.001 1.001 
1.004 1.005 1.006 1.007 1.008 1.009 
1.015 1.017 1.018 1.020 1.021 1.023 
1.038 1.040 1.043 1.046 
1.060 1.063 1.067 1.070 1.073 1.077 
1.097 1.101 1.106 1.110 1.115 1.120 
1.146 1.151 1.157 1.163 1.169 1.176 
1.211 1.218 1.226 1.234 1.242 1.251 
1.297 1.307 1.317 1.328 1.339 1.351 
1.414 1.428 1.442 1.456 1.472 1.489 
1.577 1.596 1.617 1.638 1.660 1.683 
1.815 1.845 1.876 1.909 1.943 1.979 
2.193 2.224 2.292 2.348 2.408 2.468 
2.860 2.975 3.062 3.177 3.302 3.436 
4.369 4.633 4.933 5.275 5.675 6.171 

10.87 12.94 16.00 21.13 39.31 961.9 


1.002 
1.010 
1.025 
1.048 
1.081 
1.125 
1.183 
1.259 
1.362 
1.504 
1.708 
2.017 
2.538 
3.584 
6.803 


1.003 
1.011 
1.027 
1.051 
1.085 
1.130 
1.190 
1.268 
1.375 
1.521 
1.733 
2.056 
2.611 
3.753 
7.423 


1.003 
1.012 
1.029 
1.054 
1.089 
1.135 
1.197 
1.278 
1.387 
1.539 
1.759 
2.099 
2.688 
3.937 
8.280 
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1.6 21.400 15.860 


1.7 
1.8 
1.9 
2.0 
2.1 
2.2 
2.3 
2.4 
2.5 
2.6 
2.7 
2.8 
2.9 


3.0 


3.1 


4.529 
2.386 
1.542 
1.097 
0.814 
0.625 
0.488 
0.381 
0.298 
0.233 
0.175 
0.127 
0.085 
0.048 
0.014 


4.175 
2.261 
1.483 
1.051 
0.791 
0.609 
0.477 
0.373 
C.292 
0.225 
0.170 
0.123 
0.081 
0.044 
0.010 
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TABLE 
Second Quadrant 


All values are negative 


12.540 10.350 8.800 7.646 6.732 6.022 5.432 
3.868 3.601 3.364 3.155 2.968 2.799 2.647 
2.156 2.056 1.967 1.883 1.805 1.737 1.660 
1.430 1.380 1.332 1.287 1.244 1.203 1.169 This: 
1.024 0.994 0.965 0.936 0.910 0.885 0.861 


publishe 


0.770 0.750 0.730 0.712 0.694 0.686 0.658 volume. 
0.594 0.579 0.565 0.551 0.538 0.525 0.512 
0.464 0.453 0.442 0.432 0.422 0.411 0.401 
0.363 0.355 0.347 0.338 0.330 0.322 0.314 
0.285 0.278 0.270 0.263 0.257 0.250 0.244 
0.219 0.214 0.208 0.201 0.197 0.192 0.185 


and the 


0.165 0.160 0.155 0.150 0.145 0.140 0.136 
0.119 0.115 0.110 0.105 0.101 0.097 0.093 

0.077 0.073 0.070 0.066 0.063 0.060 0.056 0.052 
0.040 0.036 0.033 0.030 0.027 0.024 0.020 
0.007 0.004 0.000 
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THEORY AND TEST RESULTS THE FATIGUE METALS 


paper was the basis for oral presentation the Joint ASCE-IABSE 
Meeting the New York Convention, October 1958. All Joint Meeting papers 
published Proceedings Civil Engineering will reprinted one 

volume. 


SUMMARY 


This paper presents empirical fatigue relationship between the mean 
stress and the alternating stress cycle; five material characteristics; 
and the resulting number cycles failure. The influence rate cycling, 
strain-hardening, and stress concentration considered separately. 


List Symbols 


upper limit stress range that leads failure after cycles— 
with stress concentration present. 


limit range nominal stress that leads failure after 
cycles—with stress concentration present. 


lower limit stress range that leads failure after cycles. 
min 


alternating stress. 


fatigue strength for complete stress reversal and cycles—with 
stress concentration present. 


Note: Discussion open until March 1960. Separate Discussions should submitted 
for the individual papers this symposium. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2222 
part the copyrighted Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. October, 1959. 


Pres., IABSE; Swiss Federal Inst. Technology, Zurich, Switzerland. 
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nominal fatigue for complete stress reversal and cycles—with 
stress concentration present. 
asymptotic value fatigue strength for complete stress reversal 
short-time tensile strength obtained conventional test. 
Caz asymptotic value creep-rupture strength for infinitely long 
time. ser 
compressive strength (Fig. 16). (a) 
relaxation stress for constant strain; for instance, stress wire (b) 
held constant length. 
initial stress material held constant strain. 
asymptotic value relaxation stress after infinitely long time (c) 
during which material held constant strain. 
invariant creep, characterizing the creep behavior material 
given temperature. 
Ther 
number identical stress cycles failure. Note from Fig. tigue: 
the moment maximum stress has been applied just prior and cra 
beginning cycling. general 
variety 
illustra 
Eq. (11) temper 
materi 
plete 
Subscripts: ure, 
for infinitely long time. 
refers nominal stress (based net area) for conditions with 
stress raiser. 
refers complete stress reversal. 
separa 


a 
a 
a 
| Ag 


hat 


INTRODUCTION 


The importance fatigue steadily increasing structural engineering 
well mechanical and aeronautical engineering. Fatigue failures cause 
losses every year, but only small portion these losses becomes 
known. great deal human effort and expenditure for tests are 

currently being devoted fatigue studies; however, the study individual 

problem represents isolated case, without general significance, unless 

can fitted into the framework comprehensive relationship. 

comprehensive relationship must furnished Theory Fatigue, which 
serves several important purposes: 


(a) provide insight into the various aspects governing fatigue behavior 
and their relationships; 

(b) permit prediction material behavior under any particular set 
cyclic stresses from the least possible number material character- 
istics, and thus permit the engineer evaluate specific problem 
without the need for special tests; 

(c) make best use available test data; 

(d) guide the researcher setting his test programs obtain the 
most significant data with minimum effort and cost. 


There are two possible avenues approach setting theory fa- 
tigue: the theoretical approach based the mechanism crack formation 
and crack propagation, and the empirical approach which attempts set 
generally valid analytical relationship between test data obtained for wide 
variety conditions and material properties. the present state the art, 
the empirical approach appears more promising. Hence, the object 
this paper present analytical relationship arrived empirically, 


illustrate examples from broad selection conditions and show that 
fits test data very closely. 


Statement Problem and Method Attack 


used this paper, the expression “fatigue failure” means failure under 


repeated applications stress that below the short-time tensile strength 


the material. (1) Generally, the stress cycle, defined the mean stress 

and the alternating stress (Fig. 1), well the rate cycling and the 
temperature, are imposed external conditions; then, the properties the 

material govern the number cycles failure. 

comprehensive theory presented here covers all ratios mean stress 

alternating stress, from the static stress-rupture case the case com- 


plete stress reversal. Furthermore, covers any number cycles fail- 


ure, from single loading (see Fig. infinity. 
Material properties change with temperature; the theory offered here ap- 
plies for any given set material properties, i.e., for any given tempera- 


ture. Hence, temperature will considered constant for any particular 


case. 
order keep the relationships reasonably simple, the rate cycling 


will dealt with separately. For the same reason, the basic theory con- 


fined metals without strain-hardening effects and without stress concen- 
tration the area failure. These influences will also considered 
separately. 
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Fig. Definition stress-cycles and stress nomenclature. 


For the conditions stipulated above, theory fatigue failure es- 
tablished provide relationship between the material properties one 
hand and the mean stress, the alternating stress, and the resulting number 
cycles failure the otner hand. The relationship reasonably simple, 
expresses the material properties minimum number readily es- 
tablished, conventional material characteristics, and fits well into test data 
covering wide variety conditions and materials. 

Fig. represents the case smooth bar steel room temperature, 
i.e., case without creep effect. The failure criteria, governed materi- 
properties, are depicted surface system coordinates: 


mean stress 
stress 


The surface divides space into region where failure should expected 
and region where failure snould not expected. The empirical relationship 
established the analytical expression this surface. 

would appear best consider first number special, simpler cases 
representing partial problems and then develop governing, over-all re- 


Partial Problems 


Complete Stress Reversal 


the case complete stress reversal (i.e., 0), the alternating 
amplitude that leads fatigue failure after cycles called the “alter- 
nating fatigue strength”, ow. The relationship between and the so- 


called Woehler curve for fatigue strength. Careful study many test results 


indicates that this relationship follows very the expression 


where 


First cycle 
On 
agty 
Fig. 
and 
sents 
for 
deper 


METAL FATIGUE 


Fig. Surface depicting the failure criteria for material with- 
out creep effect. and imposed external conditions 
re, and the material properties govern the number cycles failure. 
where Static tensile strength 
Caw asymptotic value for infinite number complete stress 
reversals 
Woehler coefficient, characteristic the material 
exponent for characteristic the material 
nship 
ses Eq. (1) becomes 
re- 
permits clear interpretation, illustrated the upper part 
Fig. ow) represents the difference between the static tensile strength 
and the alternating fatigue strength for cycles, whereas repre- 
Sents the difference between the alternating fatigue strengths for cycles and 
for infinite number cycles. The sum these differences 


i.e., constant; the relative magnitudes are governed “balance factor” that 
depends the number cycles. 
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Fig. Empirical relationship between completely reversing 
stress and number cycles failure, sential 


the other hand, Eq. (1) written the form 


j 

This expression particularly useful evaluation test results. 
relationship between and involves four material characteristics; hence, 
least four test results are required establish it. 

For the alternating fatigue strength equal the average 
short-time tensile strength and the alternating fatigue strength for 
number cycles. Since for this point and varies linearly with the 
plot vs. point-symmetrical with respect this point; hence, this 
point also the point inflection. 
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Fig. Low-cycle fatigue strength Anticorodal (Ref. 13). 


requires that i.e., initial load application prior beginning cycling 
(see Fig. for definition beginning first stress cycle). i.e., 
which has been loaded just below the tensile strength can longer 
withstand complete stress cycle with 
This conclusion can verified simple test; this test, however, re- 
careful observer and “feel” for the behavior metals just prior 
fracture. load bar just below the tensile strength, i.e., the point 
incipient necking, then relieve the load and apply the maximum compressive 
stress the bar can withstand without buckling, and finally apply tensile stress 
once more, the bar will break stress lower than that the first tensile 
loading. Fig. shows the results series such tests(12) bars 
Anticorodal these tests the ratio was -0.75 and although 
there was considerable scatter, the tensile stress failure was distinctly 
than the bar had previously withstood one several times. 


Static Stress-Rupture 


The static tensile rupture strength (Ao smooth steel bars room 
temperature essentially unaffected time; the vs. time curve es- 
straight line (Fig. 2). the other hand, steel elevated tempera- 
ture (Fig. and non-ferrous metals even room temperature (Fig. show 
decrease static tensile strength with increasing time. interest 
note from Figs. and that the shape the static stress-rupture curve 
vs. log time) very similar that the Woehler curve vs. log n). 
exactly analogous algebraic expression can set up: (4) 


(2) 


infinitely long time 

rupture coefficient, characteristic the material 


his 


Anticorodal 
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Fig. Creep-rupture strength steel 650° (W. Stauffer and 
Keller, Ref. and curve fit test data calculated empirical 


formula. 


Again, setting 


Eq. (2) becomes 


which permits interpretation similar Eq. (1b). 
Eq. (2) can also written the form 


that can linearized taking logarithms: 


log 


Fig. 

VS. 
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and Ao. 
case 
data 
applies 
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(2a) temper 
failure 


(2c) the 
value 
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Fig. Creep-rupture strength Anticorodal (Ref. 13) and 
curve fit test data calculated empirical formula. 


Fig. shows calculated from test data and the corresponding curve 


Constant Number Stress Cycles 


For constant number stress cycles failure, the variables are o,, 
Ao. empirical relationship between them has been proposed(2) for the 
case materials with creep effect. Additional study many sets test 
data has led empirical relationship for the criterion failure that also 
applies metals having creep effects: 


the invariant characterizing the creep effect the material the 
temperature under study; has the dimension stress squared. For given 
temperature, constant, irrespective the number stress cycles 


introducing the coefficients 


(4a) 
and 


the equation for the vs. curve can derived from Eq. (3) for any given 
oy: 

(5) 

particular, for the expected relationship o,, found. 
Eq. (3) can solved for for any given set values 
and that lead failure for the same number cycles: 


kg/mm? 
— 
Tests 
30 
Anticorodal 
280 
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with 
For the case Eq. (3) reduces 
from which 
setting Eq. (7), the expected relationship obtained for ma- 
terial without creep effect subjected static stress 
The two expressions for Eqs. (2a) and (7) are linked Eq. and 
the relationship 
where and Cag are the asymptotic values and respectively, for 
infinite number cycles and are obtained substituting o,y for been 
Eqs. (4a) and (4b). that 
Eq. (7) yields the stress-rupture strength for infinitely long time 
consi 
= law (Ta) it is 
Fig. shows some curves for material without creep effect 
(k2 0), and Fig. shows the analogous curves for material 
fect Both families curves are hyperbolas with common axisof 
For materials creep effect 0), all hyperbolas have 
angle -45° with the which entirely agreement with the 
already established August Woehler that “stress cycles, for which the 
difference between maximum and minimum stress does not exceed certain 
values determined experiment, even the case stresses which 
proach the ultimate tensile strength, can take place continuously, without 
rupture 
For materials with creep effect (k2 the hyperbolas intersect the 
axis points characterized The angle intersection be- 
tween the hyperbolas and the the more acute, the greater the 


q 


Fig. Empirical relationship between mean stress and alter- 
nating stress for several numbers stress cycles 
(a) for material without creep effects, 
(8) (b) for material with creep effects. 


for Incidentally, for i.e., the condition where the initial load has just 
been applied prior cycling (Fig. 1), the hyperbola becomes straight line 
expected, since for single loading failure immaterial what 
considered the mean stress and the alternating stress. this connection 


| 
(Ta) also interest note that introduction the conditions initial load- 

(3), leads 0/0. This simply means that the invariant stress charac- 

terizing the creep effect cannot obtained from ordinary tensile test. 

can shown that for the empirical law Eq. (3) the invariant 

independent the origin coordinates. This means that evaluated 

gle from test results, does not matter from where the stresses and 
were measured, providing the corresponding value introduced 
not affected). Thus, two pairs experimental values and 

mon are available, possible set 

Then, Eqs. (3), (4), (5), (6) and (7) hold for Gy, and substituted, 


respectively, for This substitution particularly 
convenient the evaluation experimental data. 
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Constant Mean Stress 


For any fixed value the mean stress expression analogous 
Eq. (1) for the stress amplitude leading failure after cycles can 
found combining Eqs. (1a), (3) and (4) into 

+ 
(10) 


where 


asymptotic stress amplitude for failure after infinite number cycles 
with mean stress 

This remarkably simple relationship between the curves for 
constant illustrated Fig. for the case Fig. 


Constant Ratio 


For any given ratio the alternating stress and the mean 
stress can expressed 


cons 
With this value for substituted, Eq. (5) becomes 
which leads quadratic equation for The 


The relationship between and the number cycles, logn, for 


& 


(10) 


les 


points 


Fig. Curves vs. for constant mean stress 


noted, however, that for this type curve the relationship Eq. (1), with 
constant exponent does not hold. 


the other hand, the value prescribed, the alter- 
nating stress becomes 


Then, Eq. (5) leads 


y a 
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The Over-All Problem 


Several combinations the special cases discussed Section may 
used describe the over-all fatigue behavior. Any such combination com- 
prises total five characteristics the material the temperature under 
consideration. 

general will found most practical and split the 
over-all problem into the problems Eqs. (1) and (3). Eq. (1) relates 
the number cycles and the material characteristics 097, and 
whereas Eq. (3) relates and for the same number cycles and 
the material characteristics and k2. The five material characteristics 
thus are k2, Caw, They may considered the unknowns 
determined test. 


Determination the Five Material Characteristics 


Tests may planned yield the five material characteristics the 
convenient and economical way. Basically, however, they can determined 
from five arbitrary combinations mean stress, alternating stress, and 
number cycles failure; for instance, from three values 
ing nj, ng, and ng, plus the two values corresponding and and 
fixed value practice, however, the evaluation should carried 
out the basis larger number tests order minimize the effect 
unavoidable scatter. Notably, the solution for the five material charac- 
teristics may based smoothed test results represented curve 
vs. and curve vs. for fixed value this case, values taker 
and for two values lead two simultaneous equations with 
and unknowns. Then, the test values vs. can used set for 
three simultaneous equations with and the remaining un- 
knowns. The problem solving five simultaneous equations for five un- 
knowns thus effectively split into solving for two and three unknowns 
time. 

The two simultaneous equations for and may solved analytically 
graphically. the graphical solution, calculated for several values 


for each equation and the intersection the two curves yields the valu 
values and k2. 


aver 

Certain data (such the following example) lends itself quite well 
analytical solution; some cases, however, the graphical solution has poss 
advantages: When more than two pairs data (for more than two values five 


are available for determining and k2, the graphical solution provides 
means for distributing errors, and thus for reducing them (triangle errors 
instead point triple intersection). Furthermore, when the available 
values for and are different for the several values ofiand 
thus not applicable, two equations the type Eq. (6) may set 
eliminate for each value The result equation between 

for each value Here, the graphical method more convenient 


the somewhat involved nature these simultaneous equations. 
Reliability the determination the five characteristics the material 

increases with the range the experimental investigation, both the 

number cycles and the spread the mean stress. This fact should and 
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Fig. Fatigue-test data for aluminum alloys 17S-T and 

24S-T (Alcoa Handbook, Ref. 14) and material characteristics 

calculated empirical formula. 

taken into account when setting test programs and when judging the relia- 
bility the calculated material characteristics. should also noted that 
set for best results more than slide rule accuracy required for many the 

analytical calculations. 

Illustrative Examples 
ally 
Aluminum Alloys 17S-T and 
Table the Alcoa Structural Handbook, 1945,(14) lists fatigue strength 
values for aluminum alloys that may regarded carefully smoothed, 

average values. These data are plotted Fig. From the smooth curves 
drawn fit the plotted points the following test data—covering the widest 
possible range for greatest reliability—were selected for the evaluation the 
ues five material characteristics: 
and 


establish and from the two pairs and values for 
8.7, first revert Eq. (9): 
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Fig. 10. Determination and 


These values now lead two simultaneous equations the type Eq. (3) 
with and unknowns: 


For 
For 8.7: 


Analytical solution the simultaneous equations leads 122.2 


sT8 


and 63.1 ksi, from which 68.1 ksi. The graphical solution 


Fig. 10. 
With determined, Eq. (6) can solved for corresponding three 
values since the respective values are known for ksi: 


Ao, 


201 150 


From and the three sets values and ow, the values and 


can calculated from Eq. (1d). Since this transcendental equation, 
must solved trial and error. most convenient make 


for calculate the corresponding three values and finally, two values 
The process then repeated until the two values are sufficiently 


close 


300 


(3) 


Fig. 11. Fatigue-test data for aluminium alloy Noral 26S-W 
(Aluminium Laboratory, Ltd., Banbury, Ref. 6), material charac- 
teristics calculated empirical formula and the corresponding 
calculated curves. 


close each other, or, other words, until the three (or more) values 
plotted against fall straight line. 
Final estimate 12.09 ksi 


Then, Eq. (4a) with substituted for ow, leads 
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Rotating beam tests 
45.6 ksi 


Tests 
°65S-W, Rep.K-RR-343- 51-30 


Fig. 12. Fatigue-test data for aluminium alloys 65S-W and 65S-T 
(Aluminium Laboratory, Ltd., Kingston, Ref. 8), material charac- 
teristics calculated empirical formula and the corresponding 
calculated curves. 
x 


from which calculated means Eq. (7a), 


=60.9 ksi 


Hence, the fatigue tests permit prediction the asymptotic stress-rupture 


strength. data are available for direct comparison; however, the re- 


duction from 68.1 ksi short-time tensile tests predicted 60.9 ksi 
for infinite time generally line with the data Fig. These data were 
obtained loading suspended bars with various weights and observing the 


time fracture.(13) 


The curves shown Fig. come very close the hyperbolas calculated 


from the material characteristics established above, can readily 
verified. 
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Fig. 13. The effect rate cycling the fatigue strength 
Anticorodal subjected complete stress reversal (Ref. 13). 


Aluminum Alloy Noral 26S-W 


Fig. shows the results obtained with wrought aluminum to- 
gether with the curves calculated the method described. Although the test 
data reveal great deal scatter, the calculated curves fit the data quite 


well; indeed, they represent natural smoothing the test results. 


Aluminum Alloys 65S-T and 65S-W 


Fig. shows data obtained from rotating beam tests (complete stress re- 
versal) two aluminum The two alloys have the same composition, 
but different heat treatment. interest note that the short-time 
tensile strength and the fatigue strength the lower range cycles con- 
siderably better for 65S-T than for 65S-W. However, the higher range 
cycles, the difference decreases and nearly vanishes. Ogw and were calcu- 
lated for the two sets data and within the accuracy the data, were found 
identical for the two alloys. Hence, full heat treatment does not increase 
the fatigue strength for infinite number complete stress reversals over 
that available from partial heat treatment. Furthermore, the function 
the same for the two alloys. 


The Influence Rate Cycling 


For metals which show creep effect, i.e., decrease rupture strength 
with time, the rate stress-cycling must obviously have influence the 
alternating fatigue strength. Fig. shows this effect for Anticorodal B.(13) 
Cycling rates were 170 cycles per minute for the slow pulsator and 6300 cy- 
cles per minute for the high-speed pulsator. the plots log 
the rate cycling appears cause merely parallel displacement the 
curve the direction the i-axis. The lower cycling rate corresponds 
somewhat smaller smaller number cycles failure for the same 
alternating fatigue strength. 
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Fig. 14. The effect strain-hardening the low-cycle fatigue 
strength Anticorodal (Ref. 13). 


The Influence Strain-Hardening Effects 


Fig. shows the results tests Anticorodal when subjected 
observed for each cycle the first few cycles; the increase more pro- 
nounced for the higher value Undoubtedly, this instance the 
fatigue behavior masked the progressive effect cold-working, i.e., 
strain-hardening. 

interest note that the case similar tests with 


-0.75 (Fig. 4), the effect strain-hardening was masked the fatigue be- 
havior. Hence, the influence strain-hardening becomes more pronounced 


Fig. shows the results fatigue tests extruded bar Aluminum 
Alloy seen that the effect strain-hardening felt about 
25,000 cycles for ratio and about 250,000 cycles 
for ratio 0.5. This would indicate that strain-hardening 


produces only temporary, unstable increase strength, rather than 
permanent and stable effect. 

The material characteristics were calculated from fatigue-test results 
the range cycles beyond the influence strain-hardening. The curves 
calculated from the material characteristics fit the test points well can value 
expected for the very considerable scatter the test data. particular 
interest the calculated curve for 1-0, which corresponds 
the creep-rupture strength. The marked decrease the creep-rupture 


strength compared with the short-time static tensile strength most natin 
ing. Evaluation numerous other tests has led the belief that this pro- 


nounced creep effect caused mainly the zinc component the alloy. 


this 
The 
tests 
results 


of 
| 
Anticorodal 


METAL FATIGUE 


89.4 ksi 


“Cold work" 


Fig. 15. The effect strain-hardening the fatigue strength 
aluminum alloy (E. Hartmann, Holt and Eaton, 
Ref. 7). 

Range Validity for Mean Stresses 

The theory fatigue presented Sections and and particularly 

Section 3c, generally applies cases where the mean stress zero 
tensile. certain extent also holds for compressive mean stress. 
However, obvious that there must limit the validity the hy- 
Fig. (or Eq. (5)) assumes greater and greater negative 
values. 
exactly analogous theory may worked out based compressive mean 
Stresses. the compression theory the alternating stress increases the 
mean stress increases algebraically, whereas for the tension theory the alter- 


nating stress increases the mean stress decreases algebraically (Fig. 16). 
The two cases reach the limit their respective validity the point where 
the alternating stress the same magnitude; generally negative for 
this point. 
The validity this approach very nicely confirmed results fatigue 
tests under axial compression conducted the University Illinois.(9) These 
results are plotted Fig. 17; for the portion governed tension, the curve 
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Governed 
Compression Tension 
Fig. 16. The ranges validity the empirical formulas based 
tests with tensile and compressive mean stress. 
calculated from Eq. (11b) with -1.0 ksi seen match the test 
very closely. 
should noted that the curve Fig. well that Fig. 
sents intersection the generalized surface Fig. with planes. Inthe 
case Fig. 16, the plane perpendicular the i-axis, whereas the case 
Fig. parallel the i-axis, but does not intersect and not 
parallel either the other two axes. 
10. The Influence Stress Concentration 
Stress raisers, such notches and holes, cause non-linear distribution 
the elastic stress over the critical section. the case single loading 
ductile materials, plastic flow occurs before failure, causing almost value 
perfect equalization stress. Hence, the tensile strength over the net area 
bar with stress raiser almost equal the tensile strength bar 
smooth bar. the case fatigue failure, the other hand, where mean 
the non-linear stress distribution plays decisive role; the nominal 
strength max (based the net area) less than the corresponding fatigue 
For any given material and number cycles, the fatigue with 
stress raiser may compared with the fatigue strength without stress raiser} 
plotting the maximum stresses max and respectively, 
vs. the mean stress abscissa (Fig. 18). cor 
Stu 
max 
VS. found very nearly straight line.(2) the case short-time 
tensile strength, i.e., for max the ratio approaches the 
value ¢=1. the case complete stress reversal, 
(13) 
The straight-line relationship may thus expressed 


(14) 
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Calculated with =750 
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366 ksi 
Fig. tests axial compression (N. Newmark, 
Mosborg, Munse and Elling, Ref. and curve 
calculated Eq. (11b). 
value cycles for steel and 106 cycles for aluminum alloys. 
addition showing the relationships between the fatigue strength 
bar with hole and smooth bar, Fig. shows the influence compressive 
mean stress and the ranges governed tensile and compressive mean stress. 
should noted that the fatigue strength welded joints follows the 
igue relationship the one shown above for elements with stress raisers.(11) 
11. Relaxation Pre-Stressed Wires 
Relaxation means decrease the stress wire that being kept 
constant length. The initial stress, “pre-stress” the final stress 
cannot exceed asymptotic value irrespective the magnitude 
Study considerable amount test data has led the conclusion that 
12) 


the relaxation process follows law analogous Eq. (1): 
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The 
Fig. 18. The influence stress concentration the fatigue 

strength steel bars subjected longitudinal force (M. and 

Eichinger, Ref. 10). 
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Fig. shows the results relaxation tests steel together 
with the curve representing the above relationship. 1951 
Relaxation pre-stressed wires particular importance the design 
pre-stressed concrete members. present the effect relaxation Nr. 
often greatly underestimated. the future, greater allowance should 
the design pre-stressed concrete members for the loss force inthe 
steel due relaxation. 
die 
Bau 


Steel 
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Fig. 19. Test data relaxation stress steel wire held 
constant length (F. Levi, Ref. 12; Wire No. 13) and curve fit 
test data calculated empirical formula. 
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published Proceedings Civil Engineering will reprinted one 
volume. 
ABSTRACT 


introduction, the penetration bending stresses due edge effects 
considered for the cases flat plates, cylindrical and rotational shells. The 
essential part the paper estimate edge effects and their penetration 
for the case hyperbolic parabolids supported straight generatrices. 


SYNOPSIS 


The structural analysis thin shells usually performed determining 
the state membrane stress due the loads and introducing bending 
corrections the neighborhood the shell boundaries. The corrections are 
required satisfy the boundary conditions which cannot satisfied 
membrane stresses alone. The penetration the bending stresses into the 
shell depends the thickness and the dimensions the shell, and also the 
characteristics the middle surface the shell. 
comparison made bending stresses due applied boundary 
moment flat plates, cylindrical shells, rotational shells, and saddle shells. 


Note: Discussion open until March 1960. Separate discussions should submitted 
for the individual papers this symposium. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2223 
part the copyrighted Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. October, 1959. 
Prof. Civ. Eng., Director, Inst. Flight Structures, Columbia Univ., 
New York, 
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The comparison illustrates the mechanism bending stress penetration 
the most commonly used types thin shells. 

The behavior the shallow hyperbolic paraboloid supported along four 
straight generatrices studied means Vlasov-type equations. The 
boundary disturbances added the membrane stresses are known 
quantitatively rather special case, and qualitatively only the general 
case. Their determination remains open problem. 


Flat Plates 


Consider simply supported, rectangular plate sides and which 
(Fig. 1). The bending moment along the axis a/2 the plate, when 
much larger than that 
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can derived means solution given Timoshenko.(1) Assuming 
Poisson’s ratio equal zero, the ratio given by: 


(1) 
where 
-mb 


Fig. gives versus y/a for should observed that be- 
comes zero 1-2/a. b/a approaches infinity, becomes first zero 
the center the plate and the penetration, defined the equation: 


(3) 
becomes 


The boundary rotation due the sinusoidally distributed boundary 
moment also sinusoidal, and for has amplitude 


Circular Cylindrical Shells 


Moment Applied Longitudinal Boundary 
Consider thin cylindrical shell radius length and thickness 


shown Fig. The values the bending moments due sinusoidal 

moment intensity applied longitudinal boun the thin cylindri- 
cal shell have been computed Parme, The values are 

given terms the parameter ah/L2 and are plotted versus s/L Fig. 


for ah/L2 0.002 and 


seen that the penetration 


(5) 


equals 0.40 and 0.90 for the two cases considered. increases with increas- 


ing values and the order 2/5 for the thinnest shells con- 
sidered practice. 


The boundary rotation was also computed Parme(5) and has 
amplitude 


(6) 
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These values are the same order magnitude that the rotation Fig. 

sidered 

Moment Applied Curved Boundary 


this case, consider semi-infinite cylindrical shell radius which 
has moment constant intensity applied free curved boundary and the 
the shell (Fig. 5). The value the bending moments the shell given 
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Fig. gives the values versus For the boundary conditions con- 


and the boundary rotation equals 


When the boundary simply supported, i.e., prevented from displacing radi- 
noted that the ratio h/a becomes smaller, the penetration and ro- 


tation are smaller different order with respect those valid for plates 
and for longitudinal cylinder boundaries. 

The maximum boundary moment and transverse shear due uniform 

outside pressure the cylinder with built-in edge have the values 


The formulas this section can applied with sufficient accuracy 
sinusoidal boundary moments sinny provided the wave length ais 


not too small compared 


Rotational Shells 


can proved(7) that the formulas for the moment penetration and 
boundary rotation due end moments obtained for the case cylinders can 
applied with satisfactory accuracy rotational shells (Fig. 7), provided: 


The opening angle not more than 20°. 

The shell thickness small enough for the penetration less than 
1/3 (or 20°). 

The shell hugs the tangent sphere the boundary region limited the 
penetration. This condition, for example, not satisfied the ellipti- 
cal shell Fig. but satisfied the toroidal shell Fig. 8b. 


For rotational shells the formulas Sec. must applied with: 


where the meridional radius the boundary and the radius right 
angles the meridional section (Fig. 7). 

order gauge the accuracy these approximate formulas may 

mentioned that(8) for hemispherical dome (‘Pp 90°) h/a 1/15 the error 
the rotation 0.50%, while for conical shell with 30° and h/a 1/160 


the error 1.23%. 


The formulas Eq. (11) for built-in edges cylindrical shells give safe 

limits for the boundary moments and shears any shell supported curved 
boundary, provided the penetration evaluated with radius equal the 

7 
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Hyperbolic Paraboloid Supported Four Straight Generatrices 


The shells the preceding section had moments applied curved 


ries. Hyperbolic paraboloids are often supported, instead, along straight 


generatrices. this case cylindrical bending theory becomes unavailable 


the simplified solution the complete stress problem the shell and the 
problem must tackled from different viewpoint. 
Let the equation the middle surface the shell defined the 


equation: 


Where the axis vertical, the and axes are parallel the sides and 
the shell, and the shell rise (Fig. 9). 

partial derivatives subscripts, let: 


{ 


the displacements the direction and the straight generatrices 
the paraboloid. The rise assumed small enough for the angle 
the generatrices considered right angle. Within the same degree 
approximation, the cosine the angle between the normal the shell and 
z-axis unity, the displacement may considered normal the shell and 
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the projection shell element the y-plane has the same area the 
The strains the and directions due and are: 


The shear strain (Fig. 10a) Vx; the shear strain ifu 
(Fig. 10b) -(2c/ab)w; hence the total shear strain given by: 


(16) 


The membrane energy due these strains expressed by:(9) 
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The bending energy the shallow shell may taken equal the energy 
plate:(10) 


Indicating the components the load per unit area plate 
surface the directions the potential energy the load is: 


Variating the total energy the system with respect and letting 
Poisson’s ratio the displacement equations for the hyperbolic 
are obtained the form: 


where: 


(21) 


Equivalent equations terms and stress function have been ob- 


tained Reissner.(11) particular solution Eqs. (20) for the case 


This solution satisfies identically the first two Eqs. (20) and reduces the 
third to: 


Eq. (23) gives ab/cEh and indicates that the shell ex- 
hibits shear deformation the type shown Fig. 11, which the external 
load does work and the deformation energy produced the boundary 
that these solutions not lead suitable membrane states stress which 
can modified local boundary corrections give practically significant 
more significant, because the only unsatisfied boundary conditions involve 
displacement and the slopes and wy. fact this solution correct 
everywhere except near the boundaries. Unfortunately not possible 
obtain the boundary disturbances for this case elementary means. The 
shell behavior clarified the following considerations. 
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Let try find complete solutions with 


and this case the third Eqs. (20) becomes: 

and shows that the shallow shell behaves like plate elastic foundation 
with foundation modulus: 

(26) 

This equation has particular (membrane) solution: 


identical with the solution mentioned above. 
order have the shell boundaries bending disturbances w(x, 
20) must added Assuming, for example, much larger than that 
wis almost everywhere function only, the equation for reduces the 
equation for beam elastic foundation: 


(28) 
with: 


bending moment the due the displacement Eq. (27) 
(22) for shell built the boundary given by:(12) 


(30) 


(23) while the penetration given by: 


(31) 
with end moment (Eq. (9)). 
ant The maximum bending stress due the moment Eq. (30) is: 
The corresponding membrane shear stress is: 
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These results confirm quantitatively the qualitative statements 


Reissner Reference the effect that: the penetration small com- 


pared b/a long h/c much less than the bending stresses 


which occur the edge zone are the same order magnitude 


membrane shear stresses. 


the solution just obtained for shell large length loaded 


substituted the three original Eqs. (20), the first and third are satisfied 


while the second violated. However, this violation can removed pro- 
viding appropriate load qy, i.e., 

(35) 

the derivative inherently small except near the boundaries 

follows that the necessary loads are importance only near the 
important notice, however, that these shears are not small. 

integrating over the width the disturbed zone one finds that, 
because 13. 


the resultant the corrective forces qab/2c. This value exactly the 
membrane shear used conventional design. The difference between the 
above solution including the forces and the actual conditions lies the fact 
that the membrane shear not applied the edge the plate should but 
acts distributed manner within narrow region near the edge. 
The complete solution the stress problem the hyperbolic 
with commensurate sides can only achieved integrating Eqs. (20) 
without the assumption that everywhere, and with physically signifi- 
cant boundary conditions. These equations seem more suitable for the 
problem than formulations involving stress function for the membrane 
stresses because some the boundary conditions involve displacements. 
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THE DESIGN FOLDED PLATES 


Eliahu 

SYNOPSIS 

Various design theories folded plate structures are critically reviewed. 
The advantages and generality the Method Particular Loadings are 
stressed and outline its application given. The method illustrated 
numerical example. 


INTRODUCTION 


Folded plates have become recent years increasingly popular; their use 
ranges from roofs, silos, staircases through such variety structures 
which the engineer’s imagination sets the only limits. The three-dimensional 
action folded plates made possible roof considerable spans with 
light and economical concrete structure, erected relatively simple form- 
work. 

therefore quite natural that the design theories folded plates have 
made during the past few decades considerable strides towards both simplicity 
and accuracy. the present paper attempt made review critically 
several theories that have been used determine the moments and stresses 
folded plates. conclusion such review the author tries point out 
the advantages and generality Yitzhaki’s Method Particular 
which appears have become far the simplest and most general way 
stress analysis for folded plates. 


The Membrane Theory 


The first design theory published(2,3) was based the assumption 
homogeneous, elastic structure which the lines intersection between the 


Note: Discussion open until March 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2229 
part the copyrighted Journal the Structural Division, Proceedings the 

American Society Civil Engineers, Vol. 85, No. October, 1959. 
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refer the theory based last assumption the membrane theory 
folded plates. The moments and stresses the structure are accordingly 
found the following procedure: 


strip unit width cut from the structure two parallel transverse 


— 


sections subjected its own dead and superimposed load. With the 
ridges assumed unyielding supports, moments and reactions the 
continuous one-way slab thus formed are then determined. This 
generally referred the “slab-action” the folded plate. 

Forces equal and opposite the reactions the slab-structure the 
ridges are then applied the “plate-structure”. These forces are re- 


solved into their components the direction the plates meeting the 


common ridge. The sum these components the plane each plate 
causing bending moment the longitudinal direction; the span 
the plate the distance between two traverses. 


The bending stresses caused the extreme fiber each plate are then 


computed. general, the longitudinal stresses the fibers two ad- 
jacent plates meeting common edge will not equal each other. 
This, course, impossible continuous homogeneous structure. 
Longitudinal shearing stresses must therefore act along every ridge 
tending equalize the stresses both fibers meeting it. follows 
from equilibrium that these shearing stresses are internally 
balanced axial forces loading each plate its edge. 


The next step determine these longitudinal forces. With each 


plate can designed for the combined moment and axial load acting 


The first design theories derived the magnitude the longitudinal axial 
forces from the condition that stresses adjacent fibers both sides any 
ridge must equal. Thus series simultaneous linear equations was es- 
tablished, their number being equal the amount ridges the folded plates 
structure. Consequently the solution the simultaneous equations was quite 
wearisome, especially for great number plates the structure. 

very substantial simplification and reduction the amount work in- 


volved finding the longitudinal forces was introduced Winter and 


Recognizing the form the simultaneous equations for the determination 


the stresses analogous those the Three Moment Theorem, Winter and 


Pei applied the Cross’ Distribution Method, which all longitudinal forces 
the ridges were directly obtained one single operation. 
Later on, this method was further simplified the introduction simi- 


lar distribution which yields directly the longitudinal stress the ridge.(5,13) 


This shortcut not only obviates the need compute the stresses from the 
moments and eccentric axial loads but also gives clear physical meaning 
the distribution, which becomes easy follow and control. additional ad- 


vantage distributing stresses rather than some imaginary forces be- 


comes obvious when certain longitudinal stresses are known the outset 
the case when applying prestressing force along any plate. Illustrations 


stress distributions under such loadings are given Yitzhaki’s 


individual plates, the ridges, not undergo any displacement. Some authors 
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The Bending Theory 


Experiments(5) well theoretical considerations prior 
have shown that the assumption unyielding ridges was unjustified most 
cases. Moreover, considering the settlement the ridges and the consequent 
rotation the plates, bending moments are caused the transverse direction 
that may radically change all stresses computed the membrane theory. 
more exact theory, taking into account the settlement the ridges, was there- 
fore developed. This generally referred the bending theory folded 
plates. 

the following the basic philosophy different approaches for consider- 

ation ridge settlements will examined. 


The Iterative Procedure 


given settlement any ridge changes the slab moments and consequently 
their reactions the ridges. But the new reactions applied the ridge will 
turn give different settlement than was originally found, i.e. further 
change these reactions. One could suggest the use iterative procedure 
which the reactions are first computed assuming unyielding supports and 
then corrected for settlements. The revised reactions applied the plate 
would then cause further change ridge settlements which turn 
requires further revision moments and reactions. 

This iterative procedure, though relatively simple use, requires not only 

substantial amount numerical work but—what more 

cases convergent very slowly not all. Moreover, can shown 
that certain cases, for small angles between adjacent plates, the iterative 
procedure actually divergent. 

iterative procedure for the determination stresses due ridge 


settlements was therefore discarded favour other approaches that treat- 
the final ridge moment the unknown. The basic difference be- 
tween various theories thus concentrated the conditions and methods from 


which the final ridge moments could derived. 


would far beyond the scope the present paper dwell the de- 
and tailed differences between various bending theories. This also not neces- 
sary here since this subject was treated quite extensively book 


Born.(8) The only thing that might stated here that the first bending 


ni- theories(9,10) developed simultaneous differential equations the fourth 
order for the determination the ridge moments settlements. The 
equations, though giving exact solution, became complex that their actu- 
office practice was extremely time consuming. 
The conditions for the determination final ridge moments were derived 
from considerations slope deflection the ridges. the folded plate 
structure with stiff connections all ridges angle change between two ad- 
plates possible. other words, the sum angle changes any 
tidge—resulting from the external loads the slabs, from the ridge moments 


and the differential settlements between adjacent ridges—must vanish. Most 
authors chose the settlement any ridge independent variable, expressed 
the transverse moment proportional its second derivative and the ex- 
ternal load function its fourth derivative. They thus obtained for each 
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ridge differential equation the fourth order involving many unknowns 
ridges the folded plates structure. 

This bending theory was considerably simplified was 
the first express the ridge deflection terms stresses and thus was 
able obtain set linear equations which yielded directly the longitudinal 
stresses well the ridge moments. Since this method relatively little 


known this country, yet was the first contain the reasoning which formed 


later the basis for more recent and simpler methods, will shortly out- 
lined the next section. 


Bending Solution 


This solution aimed determining directly one set simultaneous 
linear equations the longitudinal stresses and the ridge moments. Two con- 
ditions serve basis for the determination the unknown moments and 
stresses, namely: equilibrium and ii) zero angle change any ridge. For 
the purpose the analysis stresses due external loads are treated first 
and then unit moments are applied consecutively each ridge that capable 
resisting bending moment. 

The detailed solution carried out the following steps: 


Resolve the external load into concentrated loads the ridges (Fig. 1), 


Resolve these ridge loads into their components acting the planes 
the plates meeting the respective ridge. 
Obtain the total load acting the plane each plate, 
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Compute the plate moments, considering each plate beam resting 
two traverses: 
| M, = Pr ~ (3) 
longitudinal forces which act along the ridges, i.e.: 
n 
The moments are added algebraically, considering their sign. The index 
denotes these ridge forces the primary ones, assuming the ridge moments 
Based Navier’s assumption linear stress distribution from ridge 
ridge (Fig. 2), the longitudinal stresses are related the forces 
the equilibrium condition Eq. (5): 
where denotes the cross-sectional area plate 
Eqs. (5) can established for each ridge and solving them simul- 
taneously would yield all longitudinal stresses the structure under the 
assumption hinged ridges with zero deflection them. 


However, such assumption would not conform with the actual be- 
haviour continuous reinforced concrete structure. becomes 
therefore necessary take also into account the ridge moments and 
deflections. This done the following additional steps: 


FIG. 
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The unknown moment acts ridge and causes and the 
ridges adjacent additional reactions shown Fig. Their 
magnitude given by: 


Here 
a R = - 
n-i 
n 
(1) 
due 
Reactions cause additional plate loads AP, moments and additional 
longitudinal forces which will change the magnitude the stresses given 
Eq. (5). these stresses are considered include already the influence 
due the left side Eq. (5) will remain unchanged and only will 
the total change due will When other ridges are subse- 
quently loaded moments etc., Eq. (5) turns into: 
cally 
expr 
Here the additional longitudinal force ridge due unit 
moment acting ridge n-2. 
Every equation type (8) contains five unknown ridge moments and 
three unknown stresses the total number ridges, including the free 
edges, the number unknowns will 2m-2, i.e. stresses and m-2 
moments. Yet the number equations type (8) only therefore 
becomes necessary set m-2 additional equations. 
This easily achieved stipulating the condition that the total change 
angle any ridge, T,, must vanish: 
ples 
ink 
due 
var 
sun 
pre 
MOMENT LOADING 
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Here the change angle due the moments 


due the external loading 


the change slope due differential ridge settlements that result 
from the difference longitudinal stresses between the top and bottom the 
plate. The magnitude easily established Maxwell’s law recipro- 
cal deflecticns. The interrelation between longitudinal force and change 
slope such that the longitudinal force caused unit moment that 
point numerically equal the rotation produced unit longitudinal 
force there. Now, since the coefficients the moments Eq. (8) express 
the change longitudinal forces due unit moment, they must numeri- 
cally equal the coefficients the longitudinal stresses the equation 
expressing the rotation the plates meeting ridge (Eq. (12)): 


(12) 


thus have m-2 equations type (9) and equations type (8), 
which are sufficient for the determination all unknown stresses and 
moments The solution can tabulated and relatively simple ma- 
trix form can set for the coefficients all equations. 


Vlassow’s method and later others that were based the same princi- 
ples(11,12) provided very substantial simplification the stress determi- 
nation folded plates compared with bending theories involving simultaneous 
differential equations. However, even relatively simple structure shown 
Fig. consisting only five plates, still required the quite wearisome 
solution ten simultaneous equations. 

Another substantial disadvantage this method lies the fact that the 
ridge moments determined the solution the simultaneous equations com- 
prise once both the slab moments resulting from the external load and those 
due ridge settlements. Since the first part varies along the structure 
proportion the distribution the external load whereas the second part 
varies according the elastic line, the distribution the total ridge moments 
along the structure not known. One possibility, though extremely time con- 
set the same pattern solution described above sever- 
sections along the span the folded plate. Another possibility, normally 
preferred, assume the variation load along the structure similar 
the elastic line, namely sine wave express Fourier series. 
This, however, introduces certain error the results. 
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might pointed out, though, that the determination moments and 
stresses Vlassow’s method can not only easily tabulated, but the 
equations can also directly set matrix form with geometric data only 
coefficients. This form very easily adaptable simple solution 
electronic computer, for which this method seems especially well suited. 


Simplified Bending Theory 


With substantially the same approach shown the former section the 
determination bending moments and stresses folded plates was further 
considerably simplified.(5,13) The basic idea this approach was perform 
the stress analysis several stages, for each one which all moments and 
stresses are uniquely determined. form analogous the determination 
the redundants statically indeterminate structure the analysis 
sidesway multistory frame, primary system chosen which the 
ridges undergo deflection. unit rotation then given each plate sepa- 
rately and all stresses and deflections resulting from are computed. 
superposition, the final magnitude each plate rotation found and conse- 
quently all moments and stresses can determined. 

this way the solution reduced the determination the plate ro- 
tations only, which the structure Fig. for example, would involve only 
three simple simultaneous equations. 

Since very clear and concise description this method was given with 
numerical example recent paper Simpson, (13) there need re- 
peat here its detailed solution. 


Method Particular Loadings 


With the bending theory thus simplified might appear further 
treatment was necessary. This might have been the case were not for 
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certain lack generality inherent the methods described above which are 
not always directly applicable. 


the case, for instance, when more than two plates meet common 
ridge when there are different support conditions for individual plates—as 
for edge members supported intermediate columns—the above method deal- 
ing with plate rotations would not yield straightforward solution. Moreover, 
since the moment causing the unit plate rotation distributed along the ridge 
definite form, namely sine wave similar the elastic line, con- 
tinuous shells cannot treated. addition, take into account torsional 
stiffness edge beams shear deflection short plates and certain other 
cases would impossible the simplified theory the former section. 

Along similar lines reasoning method was developed instead 
Yitzhaki(1) that dealt directly with deflections rather than with rotations. This 


not only made possible achieve most general solution including all 
special cases, but also yielded directly the actual magnitude the vertical 


deflection along any ridge the structure. This, incidentally, many 
cases quite critical and always deserves careful checking. 

the Method Particular Loadings the total moments and stresses are 
obtained superimposing several stages loading the folded plates 
structure. The first stage deals with the external loads. clearly tabulat- 
form the slab moments, ridge reactions and plate moments are computed, 
the stresses all ridges are found and the plate deflections and slab moments 


due them are determined. These moments cause additional reactions the 


ridges, termed superfluous loads 
Additional particular loadings are then required deflect the slab 


conform with the deformed plates structure. These loadings will consist 


group vertical ridge loads that are externally equilibrium 
moment acting along ridge, both being assumed initially unity. The distri- 
bution these loads along the structure always according functions, 
called normal functions, that will yield plate deflection proportional the 
load distribution. applying the particular loading function similar 

the plate deflection the superfluous loads caused them can added direct- 
the initial ridge loads, that correction for the plate settlement due 


the additional loading required. The variation the particular loading 


along the structure can taken function similar the elastic line 
fixed beam the support conditions the folded plates are such any 
function similar the elastic line given structure. 
For each particular loading the procedure outlined above for the external 
repeated leading the determination the superfluous loads. The 


individual particular loadings are then superimposed cancel the 
primary superfluous load each ridge. This procedure again analogous 


the determination the holding forces the sidesway analysis multi- 
story frames the Cross’ Distribution Method. 
The separation the solution into the first stage dealing with the external 
load and the stages additional particular loadings permits the exact determi- 


nation all stresses along any point the structure. The variation the 
due the first stage proportional the distribution the ex- 


ternal load along the span the plates. Those due the other stages vary 

different function that practically fully expressed the first term 
Fourier series but the accuracy can increased, although entirely un- 

taking additional term. Superimposing each loading accord- 

ing the function its variation along the structure gives the exact values 

moments and stresses any point. 
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step step outline this method shown the following concise 


For 


form and further illustrated fully detailed numerical example given 


Determine the slab moments and reactions considering the structure 

one-way slab unyielding supports. 

Resolve the ridge reactions into plate loads and compute plate 
and primary stresses, assuming each plate behave 


beam. With 
Distribute these primary stresses—in method shown the 
achieve the adjusted stress every ridge, which will already point 

equal for the two extreme fibers the plates meeting the ridge. each 

Compute the ridge deflections and the fixed-end slab moments resulting there 

Distribute these moments the slab structure obtain the ridge flectior 

moments resulting from plate settlements caused the external 

Determine the ridge reactions (called the superfluous loads) due will ha’ 

moments. 
Steps are the typical procedure for each particular loading. The 
total number additional particular loadings equal the amount sign 
which the slab moments are unknown. (Such will all edges except the but 
free ones, which the moments are zero, and those next them, which The 
the moment must equal the cantilever moment the marginal plate). 
folded plate shown Fig. two additional particular loadings will require sign 
Apply the slab-structure additional particular loading “a” consist: 
ing moment acting along one ridge (ridge Fig. 4). The 
bution this moment along the structure will vary with conditions 
support, but will sine-wave for simply supported structure. The 
ticular loading R,, which the sum the initial ridge load and and 
the superfluous load. 

Repeat steps and for any further additional particular loading that 
required (mp present example). 

Superimpose the individual particular loadings, each multiplied Yit 
proportionality factor and respectively, cancel the primary Pul 
superfluous loads each ridge. This done according Eq. (13). Cr: 

denotes the particular ridge load; the first index gives the ridge and 

the second the case particular loading. 
the superfluous load ridge under the primary loading 
After determining the factors and the moments, stresses and de- 
flections any point are finally obtained superposition according 
Eq. (14) and the solution completed. 
Stress: 


q 
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ise For the determination stresses any point along the structure should 
remembered that the variation proportional the external load 
constant for uniformly distributed load) and that and pro- 


By 


portional the shape the elastic line. 


CONCLUSIONS 

Within the framework single paper course impossible show 
all applications the Method Particular Loadings. was tried though 
point out the comprehensive form the solution which can easily checked 
ige. each stage and the advantage dealing directly with the deflections. Thus 
there difficulty analyze, for instance, the folded plates structure for 

than two plates meeting common ridge; the condition equal de- 

the plates meeting the ridge will simply have satisfied. 
the marginal plate rests intermediate columns, the column reactions 
will have introduced additional particular loading unknown 

These will such cause zero deflection the marginal plate 

the columns. 


The practicing designer will find helpful not have worry about any 
sign conventions that are usually required when dealing with slope rotations, 
but able reduce the work the simplest algebraic form. 

The solution the folded plates structure the Method Particular 

Loadings thus reduced problem further complexity than the de- 
sign continuous beam. 
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APPENDIX 


Illustrative Example: The Design Folded Plate Roof the 
Method Particular Loadings 


General Data 


The folded plate roof shown Fig. A-1, that was given illustrative 
example paper Simpson, (13) will here redesigned the Method 
Particular Loadings. 

The general geometric data the individual plates that will necessary 
the course computations are first summarized Table 


Fig. 
alwa 


Moment Distribution 


Assuming strip unit width the cross section the structure act 
continuous one-way slab unyielding supports, the ridge moments are 
determined the moment distribution shown Fig. A-2. Because 
symmetry, the slab may assumed rigidly fixed ridge 

The reactions the ridges are now applied the plate-structure 
equal magnitude and opposite direction. These ridge loads well their 
resolution the direction the plates meeting the ridge are shown 
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TABLE GENERAL DATA 


ive A-3. From parts and Fig. A-3 seen that the plate loads can 
always expressed general form by: 


sary 


FIG. A-2 
MOMENTS DUE EXTERNAL LOADS 
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FIG. A-3 
RESOLUTION RIDGE 
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where the values and are already given Table The plate forces 
are tabulated column Table the Operation Table, which will 
contain all data necessary complete the design. 

With the values the moments are computed each plate by: 


and with these the initial stresses are obtained by: 


section modulus—is taken again from Table The values are 


Stress Distribution 


The free edge stresses are now distributed order determine the 
final stresses each ridge. The stress distribution shown Fig. A-4 
and the results are tabulated column Table 


*This Table the copyright Prof. Yitzhaki and reproduced here with 
his kind permission. 


OPERATION TABLE 


100 
-200 


LOADING 
ADD 
PART 


118 


October, 1959 


Determination Ridge Settlements 


The maximum displacement midspan any plate its own plane 


given 


Setting 


get 


The fixed-end slab moments due vertical ridge settlements are for 
slab whose far end simply supported, slab 


The values are given column Table 


~ 


These values are tabulated column Table II. 
The vertical displacement any ridge can seen from Fig. A-5 be: 


which becomes with the notation Table after multiplying both sides the 
equation 


FIG. A-4 
STRESS DISTRIBUTION 


Slab Moment due Ridge Settlements 


PLATE AREA 
DISTRIB. FACTOR 
CARRY OVER FACTOR 


FREE EDGE STRESS 
psi 


For sl. 
ridge 


The valu 


the 


and 


respecti 
for insta 


Thes 
ridge 
values 

The 
are obta 


—K_—— 
855 +718 +505| psi 


FOLDED PLATES 
is 2-1 b* cosy 96 * cos ¢ 


For slab fixed both ends, slab the fixed-end moment due vertical 
ridge settlement given by: 


cosy 


The values pre tabulated column Table their coefficients are 


the present example 


respectively. The fixed-end moments are given column Table 
for instance 


the 18.52 


These moments must distributed get the final slab moments each 
ridge due settlement it. The distribution shown Fig. A-6. The 
values m", the final slab moments, are given column Table 

The reactions due the moments m", resulting from ridge settlements, 
are obtained by: 


VERTICAL RIDGE DISPLACEMENT 
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FIG. 
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MOMENTS RIDGE SETTLEMENTS 


and are given volumns and Table 


The values are the superfluous loads the basic external loading; they 
will now have annulled additional particular loadings. 


Additional Particular Loading 1000 


This loading shown Fig. A-7; the variation the load along the span 
the folded plate given sin which for all practical purposes may 


considered similar the shape the elastic line the structure. 


The procedure outlined above for the external loading now repeated and 


all values are readily computed aid the Operation Table and are 


tabulated there. 
The moment for the sinusoidal loading given 
The deflection given 
that for this loading 
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The fixed-end moments due ridge settlements are calculated above 
and distributed; from their reactions the superfluous ridge loads due 
are obtained. This value added the ridge load gives the total 
particular loading due the influence mg: 
The index denotes the number ridge which this particular loading acts. 
The final values the particular loading “a” are tabulated the second 


The above outlined procedure again repeated for the particular loading 
shown Fig. A-8, and all values are tabulated the third section Table 


Superposition Particular Loadings 
they According Eq. (13) now get 


ridge 128.0 208.16a 201.12b 
The solution these equations gives: 


the span 


The final value the stresses, moments and deflections are obtained 
Eqs. (14) and are summarized the bottom section Table The ridge 
ridge for instance, which due differential ridge 

settlements thus obtained by: 


The settlement ridge becomes 


ure. 
ated and 
are 


FIG. A-& 
PARTICULAR LOADING 


aN 
rey 
{ 
TAG 
al 
where 
- 
| 


ff 


STRESSES 


PRIMARY SLAB MOMENTS 


-463 


FINAL SLAB MOMENTS 


FIG. 


are 
-826 
+707 +502 psi 
-284 
+158 
3 7 
c.SLAB MOMENTS DUE 


FOLDED PLATES 123 


Assuming 1,000.000 psi, considering long time deflection due sustained 


load, get for the vertical deflection ridge 1.85". 
final moments and stresses the center the folded plate structure 
are summarized Fig. A-9. The results are, course, identical those 


obtained Simpson’s paper. (13) 
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cells are always done remotely but maintenance and equipment set-up 


maximum maintenance and equipment set-up directly (not remotely) 
personnel might follows—although exact definitions this field are 
difficult because its newness and fluid state. 
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tective clothing for all other manipulation requirements”. 


tenant spaces; and second, the design creiteria for the main structure the 
cells reflects these basically different methods maintenance and set-up. 


Journal the 
STRUCTURAL DIVISION 


Proceedings the American Society Civil Engineers 


ALPHA-GAMMA “HOT” CELLS FOR ENTRY PERSONNEL* 


John Ruddy! 


ABSTRACT 


The basic engineering definitions and the effect the philosophy max- 
imum set-up and maintenance directly people design criteria dis- 
cussed for this type “hot” cell. example this kind cell decribed 
broad outline. 


Alpha-Gamma Cell radiation shielded and contamination 
sealed enclosed space where large quantities hazardous radioisotopes can 
handled safely. The process experimental manipulations within the 


often partly done men entering and working directly the cells after most 
the hazardous material has been remotely removed. 
reasonable definition Alpha-Gamma “Hot” Cell designed for 


“An Alpha-Gamma Cell Complex designed for operation with maximum 
maintenance and set-up directly people one which uses only proven 
maintenance-free remote handling apparatus and which uses people pro- 


Engineers responsible for the design these integrated “Hot” Cell Com- 
plexes are interested the fundamental factors affecting their safe operation 
for two main reasons. First, the degree remote versus the planned direct 
maintenance the cells affects their general arrangement and many appur- 


Note: Discussion open until March 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. Paper 2239 
part the copyrighted Journal the Structural Division, Proceedings the 
American Society Civil Engineers, Vol. 85, No. October, 1959. 

Presented the May 1959 ASCE Convention Cleveland, Ohio. 

Chf. Engr., Architectural Planning Div., Brookhaven National Lab., Upton, 
Long Island, N.Y. 
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present practice the primary design and layout responsibility for such 
“Hot” Cell Complex will most likely fall the lot engineer with either 
civil-structural with mechanical background. any event, both these 
branches engineering will substantially involved with each other well 
with electrical engineers and with nuclear specialists such health- 
physicists. These circumstances emphasize the necessity for 
and common technical language between the various branches the pro- 
fessions. 

The civil-structural engineer, because his background 
large and heavy construction concrete and steel, uniquely able design, 
specify and build this type facility. These “hot” cells, however, bear many 
resemblances large machines and because this there has real 
and reasonably complete understanding between the mechanical and civil de- 
sign engineers successful project result. 

This means, among other things, that the structural design engineer must 
able understand and sensibly apply the usual mechanical commercial 
tolerances any surfaces access holes into which mechanical 
items must mount. For example, foot shielding wall, usually 
economically impossible hold sleeves closer than 1/16" location 
bore standard construction practice and, addition, the squareness 
these holes mounting surface can deviate several degrees. Therefore, 
the mechanical alignment and fit tolerance requirements the mechanisms 
mounted inserted are more stringent than those obtainable ordi- 
nary construction, would seem advisable use rough sleeve into which 
second more accurately aligned sleeve set after the heavy work done. 
The space between the rough and finished sleeve can filled with lead wool 
some other shielding material. the other hand, the mechanical design 
engineer cannot realistically (and economically) expect machine-like ac- 
curacy heavy construction which the last analysis built contrac- 
tors—usually lump sum price basis—with laborers, iron workers and 
similar rough and ready trades people. 

There effect the type cell usuage; that is—pilot-model, re- 
search and development fundamental research—on the design 


the maximum use personnel for set-up and maintenance. The “Hot” 


Cells for fundamental research because extraordinary variations their 
uses are probably best suited the philosphy maximum personnel main- 
tenance and set-up. the operations “Hot” Cell become more repeti- 
tive over longer periods time, for example, thru the research and de- 
velopment, the pilot-model the production “Hot” Cells then does the swing 
less set-up and maintenance people become more attractive, except 
that production cells may not able economically tolerate remote hand- 
ling equipment breakdowns. The final decision how far replace- 
ing direct maintenance with remote manipulations depends the designer’s 
evaluation these three main variables: the proven dependability the 
remote handling equipment considered; the danger contamination 
personnel working the cells; and the effect the overall operation 
the failure either the above two maintenance and set-up methods. 

description recently completed “Hot” Cell for metallurgical engi- 
neering research can serve example cell complex for the acquire- 
ment fundamental knowledge and for use early research and develop- 
ment. general location the “Hot” Cell Complex for metallurgical re- 
search shown Fig. where the “Hot” Cell the semi-works area 


j pi 


either 
these 


many 
real 
vil de- 


must 
rcial 
usually 
‘ion 

chanisms 
which 
done. 
wool 
design 
ac- 
itrac- 

and 


re- 
Hot” 

their 
repeti- 
de- 
swing 

hand- 


engi- 

acquire- 

re- 
area 


ALPHA MONITOR STATION 


DECONTAMINATION 


— 


GAMMA IRRADIATION 
FACILITY 


= GRAPHY 


“ 


wast 
REACTORT wacnine SHOP 
STOCK ROOM 


u 


Hi LEVEL 
WAS'E TREATMENT 

—CcOLD 


| ume 

ce 

LECTRICAL 

| 


co.o AREA 


MECHANICAL EQUIPMENT 
ROOw 


GROUND FLOOR 


the west end the existing Hot Laboratory. large roll-up door adjoins 
the yard and provides for efficient use available space for access and 


permits handling large irradiated pieces. The cell cross-section and plan 
shown Fig. This cell designed handle radioactive material 


emitting approximately 1,000,000 roentgens per hour uranium fission pro- 


ducts its equivalent. 
There are two items about this layout which are examples important 


design criteria. The first the limitation “contamination”. This done 


subdividing the cell into smaller sections. The second the step-wise 
reduction the degree contamination. This done having separate 


| | i | 


QONTAMINATION = 
SN TRE RE 


MOD 6 MANIPULATOR 


PORT 


GAS OFF ° 
anding 
ngin 
| d MAIN FL HOT AREA 
RO || DENSITY = 250 */cu FT. 
| HIGH DENSITY 
4 LJ 
|| EWING WINDOW STEEL SHIELDING 


t 


128 October, 1959 


spaces series for the Alpha containment doors, the isolation rooms and 
the containment control area between the “Hot” Cells and the “clean” opera- 
ting area front the cell the manipulators. 

The inside dimensions the cell are feet along the operating face, 

feet deep and feet high. The front and rear walls are poured ilmenite 
aggregate concrete feet thick with density 240+ lbs per cubic foot. The 
ceiling the same material 2-1/2 feet thick. The end wall inside the 

room steel container feet out out and filled with specially graded 
dry ferrophosphorus fines with density 250+ lbs per cubic foot. This 
wall, shown Photo can removed extension the cell needed 
the future, and addition, changes the cell access can made 
after temporarily removing the dry ferrophosphorus. 

Steel work was prefabricated off-site. Unit sections were built with 
manipulator ports, window-liners and shielding plug liners place. The 
sections were assembled and welded on-site and then the walls poured the 
usual way but batches with ilmenite concrete. Some details can seen 
Photos The ceiling and roof sections were also prefabricated 
with access plug liners place. These were then erected, welded and poured 
with ilmenite concrete. After completion the cell proper, the isolation 
and shower rooms were built using non-shielding concrete blocks. The 
ceilings the isolation rooms are made removable concrete slabs the 
area directly front the Hot Cell doors. This allows the building overhead 
bridge crane used the isolation rooms. All surfaces are made de- 
contaminable keeping them reasonably smooth and coating them with vinyl 
expoxy paints. The “Hot” Cell floor large stainless steel pan, welded 
place. 

Three (3) lead glass viewing windows are installed the operating face 
the “Hot” Cell, one for each section. These are oil-filled windows with 
shielding value equal the 240 per cubic foot concrete walls. The instal- 
lation one these windows shown Photo 

The gamma ray shielding “Hot” Cell access doors are horizontally bi- 
parting and are steel inches thick and feet high. They weigh approxi- 
mately tons per set each). They hang from I-beam wheeled 
carriages and gear and chain arrangement provided for manual operation. 
They overlap the opening and effort made seal the edges. The con- 
tamination (alpha and other isotope) containment doors which seal the en- 
trance from the cell the isolation room consists two sets gasketed 
hollow metal doors with built-in windows. 

the cell provided 400 watt color improved mercury 
vapor lamps, each independently circuited. Four these lamps are recessed 
into the shielding wall above each viewing window. Each lamp fitted with 
quick disconnect fitting, requiring only one quarter turn, facilitate re- 
moval remote handling equipment. 150 watt incandescent lamps provide 
illumination when personnel are required work inside. These lamps will 
also used for emergency light case failure the mercury vapor 

| 


lights. 


The introduction ports for the entrance hot samples, filters, tools, 
etc., are modification the Argonne Tool Transfer Lock and had 
built flush with the outside wall order leave clearance for the 16" steel 
doors. disposable wiper used keep the ports both sealed and clean. 
The waste removal ports are located the isolation room side, one each 


3 
3 


129 


' 
teel 
clean. 4 
ineach 
Figure 
IS 


October, 1959 


71 ai 7 
Figure 
7 


“HOT” CELLS 131 


Figure 


room, and are standard weight 14" steel pipe mounted the cell wall. 
Shielding provided remotely operated sliding lead door 10" thick. The 
isolation room end the port recessed leave clearance for the gamma 
doors, and machined for two rings. 

shielded waste carrier designed operate with these waste ports will 
have built-in scoop push the waste can into the cell and withdraw the 
loaded can into the carrier. The waste can will have length mil poly- 
vinyl chloride tubing attached it. The other end the tubing will fas- 
tened means “O” ring the outside end the waste port. The tube 
thus sealed both ends serve contamination barrier. When the 
inside shielding door opened, the waste carrier acts shielding door 
during the withdrawal the “hot” waste can. After withdrawal the carrier 
shielding doors are closed the bag, the carrier pulled back and 
portable commercial heat sealer used cut and seal the bag. intro- 
duce new can, new bag slipped over the sealed remnant the old bag 
which removed with next load “hot” waste. effective particle- 
tight contamination seal thus maintained all times. 

pair Argonne Model manipulators will mounted above each 
the three viewing windows. Each arm will protected the latest Argonne 
designed boot which can changed without entering the cell. Four the 
six arms will have lateral rotation the slave arm relative the master 
arm. Semi-automatic control angular separation indexing standard 
all six arms. Lead shielding each arm the wall liners equivalent 
the “hot” cell wall shielding. 

auxiliary arm mounted commercial type boom will used the 
cutting section because the limited lifting capacity the Model manipu- 
lators. operated from small console “joystick” control. 

The one-ton, two-speed electric remotely operated crane bridge-type 
and will service the cutting section where the heaviest work load antici- 
pated. capable operating speeds between inches and feet per 
minute. can travel the other two sections removal the upper part 


the contamination barrier partitions which divide the cell into three 
subsections. 
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The polishing and testing sections will have commercial lift hoists in- 
stalled booms aid the Model manipulators. The booms all three 
sections are the same type order allow the hoists the auxiliary 
arm interchanged required. mock-up some the remote op- 
erated experimental equipment shown Photo 

Minor maintenance and preparations for entry into the cell will ac- 
complished remotely means the Model manipulators and auxiliary 
arm. Personnel will use one-piece plastic protective suits. This suit fea- 
tures separate face mask inside the suit with its own air-demand regulator. 
This will maintain respiratory protection for the wearer case physical 
damage the protective suit. The mask modified accept microphone 
and headset, which can plugged into the intercom system. 

Personnel, properly suited, enter the isolation rooms through gasketed 
doors the personnel change room and the door thru the partition between 
isolation rooms. The activity level the cell made low possible 
removing storing “hot” material. Modified Electro-Lux vacuum clean- 
ers remotely operated the manipulators will remove loose contamination. 
Access the “Hot” Cell proper through the gamma shielding doors and 
the set gasketed alpha barrier doors. All rooms are fitted with intercom 
jacks that personnel can maintain communication all times. 

leaving the isolation room, personnel use the contaminated shower 
(No. remove loose contamination from the suit and thorough scrub 
down done prior removal the suit from the person. After removal 
the suit, surveyed for contamination and hung the change room. 
Personnel then clean themselves the clean shower room (No. under 
health-physicist’s surveillance. The general area the “Hot” Cell Com- 
plex will under continual health-physics supervision. The “Hot” Cell and 
isolation rooms will instrumented relay general radiation level readings 


the operating face the cell which the “cold” area. 
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The ventilation system designed maintain pressure gradient, be- 
tween the “cold” operating area and the “Hot” Cells, approximately 1/2" 

water. This will accomplished placing automatic sensing devices 

the “Hot” Laboratory Building’s radioactive off-gas system. Disposable 
first stage roughing filters the cell are designed introduced through 
the introduction ports and installed the Model manipulators. Used filters 
will removed through the waste disposal port. Contamination that passes 
through these roughing filters will removed absolute filters the ex- 
isting “Hot” Laboratory radioactive off-gas system. 

The operation, more properly, the experimental manipulations that will 
“Hot” Cell described will not only manifold but more im- 
portantly, absolutely impossible complete specification even after the cells 
are built. This means that during the planning work the maximum economical 
provision must made for flexibility use. Since general human being 
one the most flexibile machines, there will always tendency 
use the maximum people rather than remote control manipulation for 
set-up and maintenance “Hot” Cells requiring flexibility. The primary 
meaning the structural design engineer these considerations that all 
the design criteria which sets prior detail design requires that 
his layout and construction should take into account both the flexibility and 
limitations both personnel and remote equipment. 

conclusion, may said that the effect structural design criteria 
for “Hot” Cells using the philosophy maximum set-up and maintenance 
people depends the design engineer’s opinion the proven reliability 
remote-handling equipment versus the contamination exposure personnel 
the cells affects their safety. general, cells using only proven 
remote-handling equipment with the operators doing all else results this 
time the simplest layouts and structure because the smaller amount 
mechanical and electrical machinery. 

Since the continued development and improvement remote-handling 
equipment unquestionably expected, obviously sets the trend. This 
trend definition will tend draw the now divergent philosophies opera- 
tion and design more closely together and will time possibly eliminate 
them. seems reasonable then say that only proven remote handling 
equipment should used “Hot” Cells used for fundamental research and 
possibly for production which means that these types would use the maximum 
people for the remainder the work which has done inside cells. 
This leaves the research-development and pilot-model cells for the field 
work the development and research required perfect the more sophisti- 
cated remote handling machinery which might well require less personnel 
set-up and maintenance the future. 


ig 


7 
J 


/2" 

ible 
ough 
filters 


will 
n- 

cells 
being 

ary 

all 
hat 
and 


teria 
rsonnel 
unt 


This 
opera- 
ate 

and 
cells. 
eld 
nnel 


October, 1959 


Journal the 


STRUCTURAL DIVISION 


Proceedings the American Society Civil Engineers 


CONTENTS 


DISCUSSION 


Design Long Reinforced Concrete Columns, Broms 

and Viest. (Proc. Paper 1694, July, 1958. Prior 

discussion: 1952. Discussion closed.) 


Generalization the Successive Approximation Method, 
Panagiotis Moliotis. (Proc. Paper 1816, October, 
1958. Prior discussion: 2043. Discussion closed.) 

Panagiotis Moliotis (closure) 


Equivalent Systems for the Deflection Variable Stiffness 
Members, Fertis and Zobel. (Proc. Paper 

1820, October, 1958. Prior discussion: 2043. Discussion 
closed.) 


Design Pier Bent and Rigid Frame Computer, 
Charles Tung. (Proc. Paper 1854, November, 1958. 
Prior discussion: 1952, 2043. Discussion closed.) 

Charles Tung (closure) 


Biaxially Loaded Reinforced Concrete Columns, Kuang- 
Han Chu and Algis Pabarcius. (Proc. Paper 1865, 
December, 1958. Prior discussion: 1988, 2075. Discus- 
sion closed.) 


Kuang-Han Chu and Algis Pabarcius (closure) 
Shear Strength Reinforced Concrete Beams, Sidney 
Guralnick. (Proc. Paper 1909, January, 1959. Prior 


discussion: 2009, 2075. Discussion closed.) 
Sidney Guralnick (closure) 


135 


141 


143 


145 


149 


Note: Paper 2237 part the copyrighted Journal the Structural Division, Proceed- 
ings the American Society Civil Engineers, Vol. 85, October, 1959. 


| 
| 
Page 
(Over) 
q 
| 
i 


October, 1959 


Strength Rivets and Bolts Tension, Munse, Petersen 
and Chesson. (Proc. Paper 1970, March, 1959. Prior 
discussion: none. Discussion closed.) 


Short Flexible Suspension Bridges for Heavy Trucks, 
Sven Olof Asplund. (Proc. Paper 2004, April, 1959. 
Prior discussion: none. Discussion open.) 

Eremin 


Analysis Two-Column Symmetrical Bents, 

Cooke. (Proc. Paper 2037, May, 1959. Prior 

discussion: none. Discussion open.) 
Glogau 


136 
Page 
159 


Page 


157 


159 


161 
161 
162 


DESIGN LONG REINFORCED CONCRETE 


Corrections Discussion J.G. MacGregor and Seiss 


Proc. Paper 1952 the February Journal the 
Structural Division, 200 line 10, please change equal 103, 
106, and 110 percent, equal 110, 106, and 103 percent, 


203, line section please change the words “the length 
effect” “the effect sustained high loads.” 


Proc. Paper 1694, July, 1958, Broms and Viest 


137 
> 
: 
is 
: ‘er 
we 
3 
| ‘ 
Lat 


e! 
al 
le 
7 
Ci 
| 
| 
7 
| 


GENERALIZATION THE SUCCESSIVE APPROXIMATION METHOD 


Closure Panagiotis Moliotis 


PANAGIOTIS ASCE.—Mr. Perng-FeiGou gives differ- 
ent way calculating distribution and carry-over factors using the slope- 
deflection equations. The calculation these factors this well known slope 
deflection equations also very interesting. 

the subject paper, calculating distribution and carry-over factors, 
identical route was followed that used the classical case calculating 
these factors framed structures with point joints and not small undeformed 
areas. 

was considered more advisable follow this route because are the 
opinion that the reader will find easier follow similar method for calcu- 
lating the distribution and carry-over factors that used the classical 
case. 

The whole study aimed calculating framed structures, whose joints are 
not points but small undeformed areas, the method successive approxi- 
mations which used designers. 


Proc. Paper 1816, October, 1958, Panagiotis Moliotis. 
Dr. Civ. Eng., Athens, Greece. 
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EQUIVALENT SYSTEMS FOR THE DEFLECTION 
VARIABLE STIFFNESS 


FERTIS and authors wish thank Mr. Perng- 
Fei Gou for his discussion their article the exact and approximate 
methods equivalent systems for the deflection variable stiffness mem- 
bers. The authors agree that the steps outlined Mr. Perng-Fei Gou, 
culminating the equivalent moment diagram, are adequate for the deter- 
mination the deflection point the system and that for this deflection, 
the steps utilizing the equivalent shear diagram and the equivalent beam, may 
omitted. However, the authors have found desirable utilize the equiva- 
lent system for the following reasons: 


Physical clarification the problem. 


Any method may used determine the deflection point, or, for 
that matter, the equation the elastic curve. 


The equivalent system greatly simplifies the approach statically 
indeterminate problems. 


moment diagram, constructed parts, utilizing the loads the 
equivalent beam adaptable for rapid calculation deflections. 


The equivalent system may used research tool. For example, 
the authors have effectively utilized the concept the equivalent system 
solve the natural frequencies vibration complex structures. The 
results this research have led greatly simplified dynamic model. 


The authors have found that, general, the error for deflection calculations 
using the approximate method equivalent systems, for minimum number 
concentrated loads, less than one percent. point interest, the 
equivalent system has negligible error when applied vibration studies. 


Proc. Paper 1820, October, 1958, Fertis and Zobel. 
Assistant Profs. Eng. Mechanics, Wayne State Univ., Detroit, Mich. 


141 
_ 
: BS 
a 
| 


} 


for 


te: 
th 


DESIGN PIER BENT AND RIGID FRAME 


Closure Charles Tung 


CHARLES TUNG,! answer the comments Messrs. 
Berg and Hoffman concerning the choice the best method computation 
for the computer, the author wishes quote his statement full from his 
original paper: 


“It generally true for any problem that the best method for manual 
computation also the best method for automatic computation. The 
reason quite simple. Automatic computation requires initial pro- 
gramming the problem, and the initial programming nothing but 

form manual computation great detail. Adopting the best method 
for manual computation will make the programming the particular 
problem considerably easier.” 


The author interested keeping the programming time any problem 
minimum. This matter simple economics and, incidentally, also mat- 
ter vital concern private engineering firm utilizing electronic compu- 
tation. long method computation can achieve the required objective, 
sophistication elegance mathematics should not influence the choice 
the method. the author’s opinion, the method using the minimum amount 

programming time the solution any problem the best method 
computation for the computer. 

Most digital computer medium scale class are basically built for fixed- 
point operations. For the solution system linear simultaneous equa- 
tions, floating-point operations are necessary. One the best methods for 
the solution these equations means matrix inversion and matrix- 
vector multiplication. general, the floating-point and matrix subroutines 
are developed and supplied the computer manufacturers. Relatively, their 
development too vast project undertaken any private engineering 
firm. the case the LGP-30, the manufacturer’s floating-point and matrix 
subroutines consist approximately 1800 commands. the author’s original 
program, still using the fixed-point operations, the solution joint moments 
achieved the use moment-distribution subroutine which consists 
approximately 300 commands. Should the method slope deflection chosen 
for the solution joint moments, the 1800-command floating-point and matrix 
subroutines would have incorporated into the program and would occupy 
major portion the memory unit the computer, thus drastically limiting 
the capability comprehensive program. Again, Mr. Hoffman pointed 
out, the slope deflection method would further require secondary solution 


Proc. Paper 1854, November, 1958, Charles Tung. 
Super. Engr., King Gavaris, New York, 
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order obtain the joint moments. true that computers can solve dif- 
ficult problems difficult ways minutes even seconds. Nevertheless, 
one does not have choose the difficult way solve problem just chal- 
lenge the computer. 

Just point information, the time when the author’s original pro- 
gram was written, the LGP-30 manufacturer had not completed the develop- 
ment the floating-point and matrix subroutines. Now, again, matter 
pure simple economics. One does not sit down and tackle the 1800- 
command subroutines, when one can find much easier way out with 300- 
command subroutine. least private consulting engineer cannot afford to. 
The author finds extremely difficult visualize Mr. Berg’s metaphor the 
method slope deflection and simultaneous equations riding over bridge, 
and method moment distribution riding ferry boat. According the 
extent programming required the two methods, Mr. Berg’s picturesque 
comparisons would seem more appropriate they were exactly reversed. 

Except for the fact that solution simultaneous equations required, 

Mr. Hoffman’s method conjugate beam definitely has its merits. Asa 
matter fact, the author has used the conjugate beam principle together with 
Newmark’s numerical procedure(1) completed program for the solution 
stiffnesses, carry-over factors, fixed-end moments, and deflections for mem- 
bers variable moment inertia. This program, mentioned the 
author’s original paper, was extension the design programs pier bent, 
rigid frame, and continuous beams. 


REFERENCE 


“Numerical Procedure for Computing Deflections, Moments and Buckling 


Loads,” Newmark, Transactions, ASCE, Vol. 108. 1943, 1161- 
1234. 
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00- Closure Kuang-Han Chu and Algis Pabarcius 

to. 

the 

the KUANG-HAN ASCE and ALGIS A.M. ASCE.—The 

sque authors would like thank Dr. and Mr. Pannell for their discussions. Mr. 

Pannell’s method very interesting and hoped that this method would 

treated more fully separate paper. 

Dr. has rightly pointed out, the main contributions the authors are 

with (i) present thorough analysis the stress-strain distribution the in- 

elastic range that could applied column sections various shapes 

and (ii) clarify the basic criteria failure such members. appears 
that precisely these points Mr. Pannell does not agree. Mr. Pannell 

bent, states: “Since the force the stress block bending and compression cannot 
defined with any accuracy and indeed difficult see that ever 
defined, for stresses are not measurable quantities then accurate de- 
termination stress can obtained from the measured strain.” true 
that stresses are not precisely measurable and the stress distribution cannot 
defined with accuracy. However, there lies the paradox, the ultimate 

kling strain (say 0.003) alone cannot used predict the failure load (including 

moment) theoretical analysis, unless the stress-strain relationship 


known, matter whether this relationship real assumed. 

The real trouble using Whitney’s stress block (as Messrs. Pannell, 
and others did) that this stress block does not represent stress strain re- 
lationship. Whitney’s block can related ultimate strain only when cer- 
tain empirical coefficients have been determined. 

illustrate the above statements, let first trace the development 
Whitney’s stress block for the case rectangular section subjected 
uniaxial bending. Assume that the stress-strain relationship trapezoidal. 
The coefficients kj, defined the ASCE-ACI report are shown 
Fig. has been shown the report that 0.5 (almost independent 
0.003, may use Whitney’s block with 0.858 (see Fig. 2). However, 
the distribution parabolic than 0.554 and 0.78. can 
seen that the values depends quite bit the stress-strain rela- 
tionship. 

Now consider the case circular section. Let try replace the 
trapezoidal stress distribution with Whitney’s stress block constant 


Proc. Paper 1865, December, 1958, Kuang-Han Chu and Algis Pabarcius. 
Associate Prof. Civ. Eng. Dept., Illinois Inst. Technology, Chicago, 
Instr., Civ. Eng. Dept., Univ. Illinois, Chicago Div., Chicago, 
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Fig. Trapezoidal Distribution Fig. Whitney's Stress 


Block 


Fig. Circular Sections 
Trapezoidal Distribution and Whitney's Block 


Centroid the 
Volume Trapezoidal 
Distribution 


the 
tic and plastic por- 
tions trapezoidal 
distribution. 


Fig.4. Section Subjected Biaxtal Bending. 
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DISCUSSION 


height 0.85 1000 and 0.003 then obviously (see Fig. 
0.858) because the part the volume the trapezoidal 
stress distribution contributed the elastic portion the distribution the 
circular section more than the corresponding part the rectangular sec- 
tion. The coefficient which indicates the location the resultant the 
compressive stress trapezoidal distribution longer equal 0.5 
nor does necessarily correspond the centroid the area occupied the 
Whitney’s block. Therefore computing the ultimate load for circular sec- 
tions, many simplifying assumptions were made Whitney and 
Hogenstad. 7)* 

similar difficulty exists using the Whitney’s block for rectangular 
sections subjected biaxial bending (see Fig. 4). For example the stress 
distribution trapezoidal, order use Au’s method with confidence, 
necessary investigate the following assumptions: (i) the trapezoidal dis- 
tribution replaced Whitney’s block with the same centroid, then 
their volume would approximately equal. (ii) The edge the equivalent 
Whitney’s block would parallel the neutral axis. (iii) The coefficient 
for uniaxial bending. Special attention should paid this assumption 
(iii) because relates the Whitney’s stress block (as represented the dis- 
tance and the ultimate strain (reached distance dg). Fortunately 
shown the authors** that with the exception assumption (ii) certain 
cases, the above assumptions are fairly good. 

would very interesting learn how Mr. Pannell has related the 
Whitney’s stress block with the ultimate strain. the author’s opinion, this 
relationship should not assumed but should established certain 
stress-strain relationships. The author’s strongly believe that order 
use the ultimate strain criteria for failure, either the stress-strain rela- 
tionship concrete compression must ascertained, the range 
variation such coefficients etc. must carefully examined. 


*Reference number given the original article, Proc. Paper 1865. 
**Discussion Au’s paper, *Ultimate Strength Design Rectangular 
Concrete Members Subjected Unsymmetrical Bending” Kuang-Han Chu 


and Pabarcius, Journal the American Concrete Institute, Sept. 1958 
1253. 
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SHEAR STRENGTH REINFORCED CONCRETE 


Closure Sidney Guralnick 


SIDNEY GURALNICK,! A.M. ASCE.—The author would like take this 
opportunity thank Messrs. Chu, Warner, Cowan, Goschy and Balazas for 
the time and effort they have expended writing discussions his paper. 
The comments made these discussions certainly form significant contri- 
bution the value the paper. 

Dr. Chu has raised three questions considerable importance and at- 
tempt answer them follows. 

The first question asked was: “Why used predict the occurrence 
the first diagonal tension crack?” answer this, although failure 
“pure shear” plain concrete almost always observed “diagonal- 
tension” rupture, correct use the maximum shear stress, which 
corresponds the failure state stress describe same (since sim- 
ply the radius the pure shear stress circle failure). has been stated 
Bijlaard(3) and others, the radius the “pure shear” stress circle must 
equal magnitude the diameter the “pure tension” stress circle 
plain concrete. However, straight-line approximation the Mohr en- 
velope adopted, the magnitude the maximum shear stress failure for 
element subjected pure shear state stress given 
geometrically-derived relationship which the author’s Eq. 13. The magni- 
tude given this equation always less than the magnitude the pure 
tensile strength, and Eq. therefore represents conservative estimate 
the pure shear strength plain concrete (as was demonstrated SectionIV 
the paper.) 

The second question concerns the author’s use average stresses. The 
average normal and shear stresses were used the author’s derivations 
simply because impossible the present state knowledge predict 
the “exact” stresses the compression zone the beam after flexural ten- 
sion and diagonal-tension cracking have occurred. Consequently, recourse 
was had the average stresses even though was recognized that Mohr’s 
circle properly refers the state stress point. Since most cases 
the normal and shear stress distributions are convex outward, there 
probably not great deal difference between the average and the maximum 
stress. Furthermore, quite likely that the maximum normal stress will 
occur level other than that which the maximum shear stress occurs. 
Therefore, the location within the compression zone which the combination 


Proc. Paper 1909, January, 1959, Sidney Guralnick. 
Assistant Prof., Civil Engr. Dept., Inst. Technology, Chicago, 
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normal stress and shear stress such cause failure will very 
probably one which either one both the stress components are 
lower than their respective maximum values. Thus the inexactitude inherent 
using the average stresses probably quite low. 

The third question raised concerned the applicability the author’s 
theories the case beam section subjected partly elastic and partly 
inelastic stress distribution. Or, other words, the theory applicable 
the maximum compressive stress failure expected somewhere 
the range between fifty and one hundred percent The author believes 
that his theory applicable such cases one can correctly compute the 
quantity was pointed out Mr. Warner, necessary use 
compatibility condition obtain analytical expression for ky. Unfor- 
tunately, explicit compatibility condition used predicting shear 
strength yet definitively established. and when such compatibility 
relation established, can easily incorporated into the author’s theory. 
view the present lack explicit compatibility relation, one must 
somehow estimate when expected that the compressive stress 
failure will lie the range between fifty and one hundred percent One 


way estimate for this case assume that tne stress the tension 


steel the yield point and given the well-known equation from 
flexural strength analysis 


which and are given Eq. 19. The minimum, and therefore the 
most conservative value for will result one uses this procedure. 

reviewing the “restrained” beams given Table-1, the author used 
slightly different approach computing ky. Since large amounts com- 
pression reinforcement were present all the restrained beams, the most 


conservative estimate resulted one used the following equation 
compute 


which the yield strain the tensile reinforcement and the ulti- 
mate compression strain for plain concrete. The value was obtained 
from the equation given Hognestad, Hanson, and McHenry(11) as: 


0.004 


The author fully agrees with Mr. Warner that further theoretical and ex- 
perimental studies are needed regard the problem formulating 
rational compatibility condition for shear failure. Furthermore acutely 
aware the pressing need for studies the behavior beams loaded 
through side stubs since these more closely represent the conditions 
found actual cast-in-situ structures than the typical top-surface-loaded 
laboratory test specimen. 

The constructive evaluation the author’s work Professor Cowan 
most appreciated. very gratifying that Professor Cowan finds many 
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DISCUSSION 


areas agreement with the author’s arguments. The one area apparent 
disagreement concerning failure criteria more one semantics than one 
physical principle. Professor Cowan argues that the tensile and com- 
pressive strength concrete are two separate physical constants the ma- 
terial, then two separate failure criteria must used. the author’s opin- 
ion, this just another way stating that least two physical constants are 
necessary completely specify the envelope circles representing failure 
states stress (i.e. the “rupture envelope” ‘intrinsic curve”) and, the 
case plain concrete, these physical constants are the tensile and the com- 
pressive strengths. Professor Cowan’s dual failure criterion differs from 
the author’s failure criterion only the details concerning the method 
specifying the exact shape the rupture envelope. Results obtained from 
either approach should comparable. 

unfortunate that Professor Cowan misinterpreted the author’s Fig. 
thus leading him conclude there was great physical difference between 
his and the author’s failure criteria. Fig. was deliberately not drawn 
scale since its sole purpose was illustrate the various geometrical quanti- 
ties and relations used the development Eq. 14a. single fixed value 
the angle was used the author’s work (as erroneously assumed 
Professor Cowan.) The angle was treated perfectly general quantity 
whose functional relations are given Eq. 11. Perhaps more meaningful 
relation for the angle results one transforms the second Eqs. 


From this equation obvious that the ratio tensile compressive 
strength, decreases angle increases. Thus seen that, contrary 
Professor Cowan’s implications, significant desparity exists between the 
author’s work and that Balmer(2) Richart, Brandtzaeg and Brown.(18) 
The author apologizes for the typographical error which appears page 10, 
reference (54) should reference (37). 

Many thanks are due Messers Goschy and Balazas for their independent 
verification the author’s Eq. 14a which central his entire theory 
shear strength. Their use the concept “reduced stress” very interest- 
ing and may have important applications building code provisions and design 
procedures the author’s theories gain wider acceptance. believed that 
the questions raised their remarks are fully answered elsewhere this 
discussion. 

view the questions raised concerning certain the assumptions made 
the author and check the validity his theory, further analysis 
presented the shear strength beams reinforced with vertical stirrups. 

Eq. 36, which based empirical expression developed Laupa, Siess 
and gives the total shear resistance, beam reinforced 
with vertical stirrups. analytic expression for which very similar 
form Eq. may derived follows. 

Considering the free-body diagram portion beam reinforced with 
vertical stirrups shown Fig. 16, the following assumptions may made: 


For the purposes calculating ultimate strength, only the stirrups 
crossed diagonal-tension crack are effective contributing the 
shear resistance the beam. 
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failure the stress the stirrups crossed diagonal-tension 
crack equal the yield stress the steel. 


the free-body diagram shown Fig. equilibrium, then 


FIGURE 


which the contribution the shear resistance the concrete and 
that the stirrups. may also written 
the stirrups the region crossed diagonal-tension crack are equally the 
spaced, then 
which 
Therefore, 
Substituting given Eq. into Eq. equation for obtained 
or, 


(41) 


cot 
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(38) 


(39) 


(40) 


(41) 


DISC USSION 


the right-hand side Eq. 41a multiplied the ratio which 
equal unity, the equation for becomes 


the angle less than 45° (the usual case*), cot greater than unity. 
Therefore conservative estimate results set equal unity 
Eq. 42. this done Eq. becomes 


(43) 


Substitution given Eq. into Eq. yields equation for the 
shear resistance reinforced concrete beam 


Eq. seen remarkably similar the empirically-derived Eq. 


One obvious advantage Eq. over Eq. that dimensionally cor- 
rect. 


Eq. may further simplified substituting for the ratio ky)/ky 
the quantity this done Eq. becomes 


the usual case less than 0.5 and therefore conservative value for 
unity. those unusual cases which exceeds 0.5, Eq. must 


used compute shear resistance. Or, other words, Eq. may always 
used provided that 


or, 


whichever the smaller. 


*It possible for angle greater than 45° but such will the case 
and only the M/Vd ratio less than about 2.5. 
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demonstrate the validity Eq. 45, ultimate strengths the restrained 
beams reported Table are compared with values predicted shear 
strength. The results this comparison test and theory are reported 
Table-4. addition, comparison test and theory for four new beam 
whose properties are given Table-3, are also reported 
Table-4. Beams and IID-IM had plain webs while the other ten beams 
reported Table-4 had web reinforcement consisting equallyspaced 


vertical stirrups. 


* Flange Width of Tee-Beam 


** Stem Width of Tee-Beam 


Test Results 


Table - 3 


Properties of Series-M Beam Specimens 


(psi) 


(psi)| (ksi) 


Table - 4 


Comparison of Beam Test Results with Theory 


Results of Calculation 


*These specimens are designated series-M beams. 


2- 


$5 a3 1/2 
2 - $7, 2 - fio 
2-07,2- #4 83 -1/2 


0! 
Beam L be piece da M/Va Long. Wed 
seam Type 
Wo. Failure Ye Yeot Pred. Error = 
x 100 
(kips) (kips) (kips) (kips) 


e 
nN ~_ w 


DISCUSSION 155 


may observed Table that the theory correctly predicts the mode 
failure all but one case. The exception beam which had low 
computed flexural strength, Beams and which failed 
flexure also had low computed flexural strengths. This underestimation 
flexural strength was due the fact that the tensile reinforcing steel used* 
had little “yield plateau” and entered the strain hardening range almost 
immediately upon reaching its yield stress. The flexural strength calcula- 
tions were, course, simply based attaining the yield stress the tensile 
steel failure and they did not reflect any extra strength which might occur 
tensile steel stresses progressed into the strain hardening range. 

Considering all beams other than the three mentioned above, the extreme 
right-hand column Table-4 shows that the theory conservative all but 
one case. Furthermore, six the beams which failed shear the percent 
error between test and theory confined the range 7.7% 15.8%. the 
remaining three beams the percent error rather high, approximately 24% 
two beams and 41% the third beam. Since the theory errs the con- 
servative side these cases and since these results may have been affected 
errors experimental technique, believed that these last three beams 
not cast any serious doubt the general validity the author’s 
theories. 

conclusion the author feels that promising analytical approach the 
problem predicting the shear strength reinforced concrete beams has 
been initiated. Much further analytical and experimental work remains 
done; however, believed that the very least the theories developed 
herein comprise rational point departure for future research. 


*All beams were reinforced with high-strength alloy steel deformed bars 
having nominal yield point 80,000 psi. 
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STRENGTH RIVETS AND BOLTS 


gratulated another fine report long series dealing with the study 
fasteners. stated that some design specifications still prohibit the use 
structural connections which are subjected direct tensile loads. 
understand that there are now only few structural specifications which 
not sanction the use rivets and bolts, either, connections subject 
tension. Bolts tension probably are used greater quantities struc- 
tures, such multi-story buildings, than anywhere else. Multi-story 
buildings, and others, have used rivets tension before the event the high- 
strength bolt, and are still doing so. 

There are many applications bolts which are used tension, for ex- 
ample the flange bolts studs high-pressure and high-temperature piping 
systems; cylinder-head studs engines and many others use for many 
years. 

stated that the cross-sectional area driven diameter rivet 
about 0.518 sq. in. compared the stress area bolt being 0.334 
in., and that lower strength rivets may develop load-carrying capacity 
approaching that high-strength bolts. This rivet area based 
rivet completely filling hole when hot driven. This condition only 
true cases very short grip lengths, and even then dependent the 
human element driving any particular rivet. Design would require the use 
the area the nominal diameter the rivet, 0.442 sq. in. for 3/4" 
rivet. 

noted that nut stripping failure occurred some the tests and 
feel should pointed out that the American Society Testing Materials 
took steps 1957 increase the nut proof load and that the current A.S.T.M. 
Spec. A325-58T now requires the minimum nut stripping values approxi- 
mately 25% higher than those the minimum bolt strength values. 


Proc. Paper 1970, March, 1959, Munse, Petersen and Chesson. 
Bethlehem Steel Co., Bethlehem, Pa. 
Bethlehem Steel Co., Bethlehem, Pa. 
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SHORT FLEXIBLE SUSPENSION BRIDGES FOR HEAVY 


EREMIN, ASCE.—The author has stated that heavy truck loading 
produces large, around inches, deflection. However, Mr. Asplund stated 
that the trucks are easily overcoming the resulting heavy grades and also that 
the truck drivers not notice the heavy grades. Yet, Mr. Asplund then 


stated that “the translations and rotations deck structure need considerable 
attention.” 


would interesting know more about the cost continuous inspection 
and maintenance the deck structure, details their original condition, and 
the cost regulating the traffic across the bridge. 

For the cast short span suspension bridges considerable study was made 
methods reducing the live load deflections. Figs. (a) and (b) are 


Proc. Paper 2004, April, 1959, Sven Olof Asplund. 


Assoc. Bridge Engr., Bridge Dept., California State Highways, Sacramento, 
Calif. 
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shown short span suspension bridges which the moving live load deflections 
are reduced hinges and tie rods. Dead load hinge Figs. (a) and (b) 
may eliminate the negative live load deflection that section. Furthermore, 

the height towers the suspension bridges Figs. (a) and (b) may 
also reduced minimum. 
The economic value short span suspension bridges important. There- 
fore, the author has made valuable contribution showing details and design 
method typical suspension bridge used Swedish hydro-electric plants. 
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ANALYSIS TWO-COLUMN SYMMETRICAL BENTS 


and Michael Jordan 


Cooke has presented interesting method which 


has been found very useful practice. does not appear aware 
however that almost identical method was published Naylor the 
“Structural Engineer” London, April, 1950. 

addition the cases dealt with Mr. Cook, Naylor applied the method 
Vierendeel girder with similar top and bottom chords. 

Further papers this subject have appeared German publications, the 
most interesting being one Ehlers and published 
Stahlbeton” Jan.1957. Ehlers shows how multibay, multistorey frame can 
dealt with equivalent two column symmetrical bent provided the 
column stiffnesses are certain ratio, which function the beam 
stiffnesses. This ratio (but not necessarily the sizes members) must 
constant for all storeys. practice particularly with concrete buildings this 
ratio can easily obtained. 

There are course large number building frames where approxi- 
mate analysis leads sufficiently accurate results when seen the light 
uncertainties lateral loads and other factors such uneven settlement. 


PEI,2 ASCE.—The author presented some numerical examples 
two-column symmetrical bents analysed the Cantilever Moment Distribu- 
tion Method Grinter and Tsao.(1) There appears nothing essentially 
new the present paper. That the Cantilever Moment Distribution Method 
applicable only symmetrical bents subjected anti-symmetrical loadings 
was pointed out the writer.(2) The extension the method un- 
symmetrical loadings first separating the given loading symmetrical 
part, together with numerical example was also pointed out the writer. (2) 

When applicable, the Cantilever Moment Distribution Method the fastest 
method known. The popularity the method will undoubtedly increase with 
time becomes better known the profession large. order avoid 
its misuse, important state the following limitations the Cantilever 
Moment Distribution Method: 


The structure must symmetrical. 
The loading must anti-symmetrical. 
The structure must two-column bent, two-chord Vierendeel 
truss. other words, section across the structure, taken normal 
Proc. Paper 2037, May, 1959, Cooke. 
Associate Prof., The City College New York, New York, 
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the axis symmetry, must not cut more than two members the 
structure. 


The last limitation inherent the derivation the Cantilever Moment 


Distribution Method. far know, has never been clearly stated any 
published paper before. 


REFERENCES 


“Joint Translation Cantilever Moment Distribution,” Grinter 
and Tsao, Transactions, ASCE, Vol. 119, 1145, 1954. 


Discussion Translation Cantilever Moment Distribution,” 
Pei, Transactions, ASCE Vol. 119, 1954, 1208. 
MICHAEL JORDAN, A.M. ASCE.—One the examples the paper treats 
the problem applying the Cantilever Moment Distribution method bents 
having variable moments inertia. The author determines the modified stiff- 
ness equating the M/EI diagram equal unity. Where the members are 
haunched otherwise tapered the expression for the area under the M/EI 
diagram will become more involved. 
method for determining the modified stiffness any member can de- 
veloped using only the ordinary Moment Distribution constants. Let 
equal the modified stiffness member AB. 
The following relationship well known. 


similar relationship exists for the columns. 


Designer, O’Neil Co., Oakland, Calif. 


any 
So: 
The coefficient can also used determine fixed end moments for 
given Story Moment, 
tiff- 


| 
ik 


Ti 
| 
* 
4 
) 


— 


PROCEEDINGS PAPERS 


The technical papers published the past year are identified number below. Technical-division 
indicated abbreviation the end each Paper Number, the symbols referring to: Air 
Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway (HW), Hy- 
draulics (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (SA), Soil 
Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterways and Harbors 
(WW), divisions. Papers sponsored the Department Conditions Practice are identified the sym- 
(PP). For titles and order coupons, refer the appropriate issue “Civil Beginning 
with Volume (January 1956) papers were published Journals the various Technical Divisions. 
locate papers the Journals, the symbols after the paper number are followed numeral designating 
the issue particular Journal which the paper appeared. For example, Paper 1859 identified 
1859 (HY7) which indicates that the paper contained the seventh issue the Journal the Hydraulics 
Division during 1958. 

VOLUME (1958) 


OCTOBER: 1790(EM4), 1791(EM4), 1793(EM4), 1794(EM4), 1795(HW3), 1797(HW3), 
1798(HW3), 1801(HW3), 1802(HW3), 1803(HW3), 1804(HW3), 1805(HW3), 1806(HY5), 
1816(ST6), 1817(ST6), 1820(ST6), 1821(ST6), 1822(EM4), 1823(PO5), 1824(SM4), 


NOVEMBER: 1833(HY6), 1834(HY6), 1835(SA6), 1836(ST7), 1837(ST7), 1838(ST7), 1839(ST7), 1840(ST7), 
1841(ST7), 1842(SU3), 1843(SU3), 1844(SU3), 1845(SU3), 1846(SU3), 1848(SA6), 1850 


(ST8), 1868(PP1), 1869(PP1), 1870(PP1), 1871(PP1), 1872(PP1), 1873(WW5), 1876 
1886(PO6), 1887(PO6), 1888(PO6), 1889(PO6), 1890(HY7), 1891(PP1). 


VOLUME (1959) 


JANUARY: 1892(AT1), 1893(AT1), 1895(EM1), 1896(EM1), 1897(EM1), 1898(EM1), 1899(HW1), 
1900(HW1), 1901(HY1), 1902(HY1), 1903(HY1), 1904(HY1), 1906(PL1), 1907(PL1), 
1909(ST1), 1911(ST1), 1912(ST1), 1913(ST1), 1914(ST1), 1915(ST1), 
1927(HW), 1928(HW1), 1929(SA1), 1932(SA1). 


1941(ST2), 1942(ST2), 1943(ST2), 1944(ST2), 1945(HY2), 1946(PO1), 


1959(CO1). 


MARCH: 1960(HY3), 1961(HY3), 1962(HY3), 1964(IR1), 1967(SA2), 1968(SA2), 
1969(ST3), 1970(ST3), 1971(ST3), 1972(ST3), 1973(ST3), 1974(ST3), 1975(ST3), 1976(WW1), 


1979(WW1), 1980(WW1), 1981(WW1), 1982(WW1), 1983(WW1), 1984(SA2), 1986 


APRIL: 1990(EM2), 1991(EM2), 1992(EM2), 1993(HW2), 1994(HY4), 1995(HY4), 1996(HY4), 1997(HY4), 1998 
(SM2), 1999(SM2), 2000(SM2), 2001(SM2), 2002(ST4), 2003(ST4), 2005(ST4), 2007 


2014(AT2), 2015(AT2), 2016(AT2), 2017(HY5), 2018(HY5), 2019(HY5), 2020(HY5), 2021(HY5), 
2023(PL2), 2024(PL2), 2025(PL2), 2026(PP1), 2027(PP1), 2028(PP1), 2029(PP1), 2030(SA3), 2031(SA3), 
2032(SA3), 2033(SA3), 2034(ST5), 2035(ST5), 2038(PL2), 2039(PL2), 


JUNE: 2048(CP1), 2049(CP1), 2050(CP1), 2051(CP1), 2052(CP1), 2053(CP1), 2054(CP1), 2055(CP1), 2056 
(HY6), 2057(HY6), 2058(HY6), 2060(IR2), 2061(PO3), 2062(SM3), 2063(SM3), 2064(SM3), 2065 


(ST6), 2066(WW2), 2067(WW2), 2068(WW2), 2069(WW2), 


(SA4), 2088(SA4), 2089(SA4), 2090(SA4), 2091(EM3), 2092(EM3), 2093(EM3), 2094(EM3), 2095(EM3), 2096 
2105(AT3), 2106(AT3), 2107(AT3), 2108(AT3), 2110(AT3), 2111(AT3), 


2114(AT3), 2115(AT3), 2116(AT3), 2118(AT3), 2119(AT3), 2120(AT3), 
2123(AT3), 2124(AT3), 


AUGUST: 2127(HY8), 2128(HY8), 2129(HY8), 2131(PO4), 2134 


2186(ST7)°, 2187(PP2), 2188(PP2). 


OCTOBER: 2190(AT4), 2191(AT4), 2192(AT4), 2193(AT4), 2194(EM4), 2195(EM4), 2196(EM4), 


2197(EM4), 2198(EM4), 2199(EM4), 2200(HY10), 2201(HY10), 2202(HY10), 2203(PL3), 2204(PL3), 2205 
(SM5), 2216(SM5), 2219(ST8), 2220(EM4), 2222(ST8), 2223 
(ST8), 2224(HY10), 2225(HY10), 2229(ST8), 2230(EM4), 2231(EM4), 


Discussion several papers, grouped divisions. 
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